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ABSTRACT 

A  study  was  made  of  Dual  Input  Describing  Functions  (DIDF)  for  non¬ 
linear  elements  with  a  view  toward  the  synthesis  problem  where  the  charac¬ 
teristics  of  the  DIDF  are  specified  a  priori.  The  study  included  a  literature 
survey  and  an  analytical  investigation  of  the  DIDF. 

Improved  methods  for  calculating  DIDF's  were  sought.  The  problem  of 
defining  the  DIDF  in  such  a  way  that  it  is  valid  for  multivalued  nonlinear 
elements  was  also  considered  and  one  method  of  solution  is  proposed.  The 
effect  of  changes  in  the  secondary  signal  waveform  on  the  DIDF  for  nonlinear 
elements  was  also  investigated. 
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CHAPTER  I 

INTRODUCTION 

s 

I 

i 

t 

t 

A  typical  nonlinear  element  N,  as  encountered  in  automatic  control 

t 

I 

applications,  may  be  characterized  mathematically  as  a  nonlinear  operator 
which  acts  on  a  scalar  input  signal  e  and  produces  a  scalar  output  signal 
y  =  N(e)  where,  in  general,  N(e)  is  a  nonlinear  function.  Iomanycases,  N(e)  is 
a  multiple  or  even  infinitely  valued  "function"  and  may  possess  a  numoer  of 
simple  jump  discontinuities.  Hereafter  the  function  N(e)  is  referred  to  as  the 
"characteristic"  of  the  nonlinear  element  N. 

One  method  commonly  used  to  analyze  electrical  networks  and  feedback 
control  systems  containing  such  nonlinear  elements  is  the  Method  of  Describing 
Functions  (DF) .  This  method  consists  of  a  linearizing  process  whereby  the 
nonlinear  operator  N(e)  is  replaced  by  a  •;  possibly  complex)  parameter 
dependent  linear  operator  called  the  describing  function.  The  DF  method 
originated  in  the  sinusoidal  analyses  of  feedback  control  systems  containing 
nonlinear  elements  and  was  therefore  originally  developed  only  for  sinusoidal 
inputs.  The  DF  for  that  case  can  be  explained  by  Figure  1. 1.  The  constant 
(possibly  complex)  gain  K(A)  is  chosen  so  that  the  output  Aj  sin  (cut  +  </>])  of 
the  linearized  representation  is  precisely  equal  to  the  fundamental  component 

« 


‘MM**,  k*-**’' 
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of  the  actual  output  y(t)  of  the  nonlinear  element,  the  latter  being  determined 
by  an  ordinary  Fourier  Series  analysis  of  y (t) .  The  remainder, 

R(t)  =  A2  sin  2o>t  +  B2  cos  2o>t  +  —  ,  of  the  actual  output  y(t)  is,  in  effect, 
neglected. 


e  =  A  sin  cot 

N(e) 

y*TJ  -  A j  Sin  COT  T  D j  COS  COT 

+  s‘n  2cot  +  E&2  cos  2cot  +  a^  sin  3cot  +  .  .  .  . 

(a)  Actual  Output  of  Nonlinear  Element  N 


A  sin  cot  | 

1C(A) 

A^  sin  (cot  +  f  j)  -  A^  n  cot  +  Bj  cos  cot 

(b)  Linearized  Representation  of  Nonlinear  Element  N 
Figure  1.1.  Replacement  of  the  Original  Nonlinear  Element  N  by  Its  DF  K(A) 


Thus,  the  approximation  of  N  by  its  DF  is  useful  primarily  in  applica¬ 
tions  where  the  signal  y < t)  subsequently  passes  through  a  filtering  process 
such  that  the  contribution  of  R(t)  at  the  filter  output  is  negligible.  Ir  fact,  it 
was  situations  of  this  type  which  prompted  the  original  applications  of  the  DF. 
As  the  filtering  action  more  closely  approximates  that  of  a  perfect  low  pass 
filter  [low  pass  with  respect  to  the  fundamental  frequency  oj  of  y(t)  ] ,  the  DF 
approximation  of  N  becomes  more  exact.  It  is  remarked  that  the  ordinary  DF 
analysis  is  valid  only  if  the  nonlinear  element  output  y ( t)  has  zero  average 
value  and  the  fundamental  component  of  y(t)  has  the  same  frequency  u>  as 


2 


the  input  e(t)  =  A  sin  cut. 


The  intuitive  (original)  approach  to  the  DF  analysis  described  above 

A 

leads  one  to  the  following  definition  of  the  parameter  dependent  gain  K(A) 

K(A)  =  1  =  Kr  +  jKj  ,  (1.1) 

^  _ 

where  K(A)  is  defined  in  complex  variable  notation  (j  =  v-l) .  The  terms  Aj 

and  Bj  are  the  magnitudes  of  the  fundamental  in-phase  and  quadrature  compo¬ 
nents  of  the  output  y(t) . 

The  DF  described  above  can  also  be  derived  by  means  of  a  seemingly 
different  approach  to  the  problem.  In  this  approach  the  unwanted  harmonics 
R(t)  are  minimized  in  the  root-mean-squared  (RMS)  sense.  In  other  words, 
one  seeks  the  function  K(A)  which  gives  the  minimum  RMS  value  of  R(t)  in 
Figure  1.  2. 


It  is  not  difficult  to  show  that  the  minimum  value  of  R  results  when 

1  2n 

K  =  — —  f  N  (A  sin  cot)  sin  cot  dcot  (1.3) 

R  ttA  J 
o 

-  i  f 

K,  =  —  |  N(A  sin  cot) cos  cot  dcot  .  (1.4) 

I  7rA  “ 

0 

The  result  [Equations  (1.3)  and  (1.4)  ]  gives  precisely  the  same  DF 
previously  described  and  illustrated  in  Figure  1.1.  However,  the  derivation 
of  DF's  by  minimization  of  the  RMS  value  of  R(t)  provides  a  convenient  means 
for  accommodating  more  general  (non-sinusoidal)  periodic  inputs  to  N. 

The  conventional  Dual  Input  Describing  Function  (DIDF)  is  defined  for 
an  input  e  which  is  the  sum  of  two  (possibly  independent)  components  e.  where 
the  e^  are  typically  chosen  as  constant,  sinusoidal,  or  random  noise  signals. 

In  this  case,  a  DIDF  or  "equivalent  linear  gain"  can  be  defined  for  each  input 
component  e^(t) .  Figure  1. 3  illustrates  how  a  nonlinear  element  with  two  sine 
wave  input  components  can  be  replaced  by  two  DIDF's. 

It  is  useful  at  this  point  to  introduce,  in  a  simplified  way,  the  underlying 
principle  of  the  DIDF  for  a  nonlinear  element  with  two  independent  input  compo¬ 
nents.  A  more  complete  mathematical  description  of  the  DIDF  and  one  means 
of  deriving  it  for  deterministic  input  components  is  given  in  Chapter  II. 

If  the  primary  component  of  the  input  to  the  nonlinear  element  is  denoted 
by  ej  and  the  secondary  input  component  is  e2,  the  DIDF  for  the  component  ej 
is  defined  as  the  coefficient  of  the  term  ej  in  the  output  y(t) ,  where  y(t)  is 
written  in  the  form 
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y(t)  =  K,ei(t)  +  R(t) 


(1.5) 


Figure  1.3.  Example  Showing  Nonlinear  Element  Replaced 

by  Its  DIDF 

and  R(t)  does  not  contain  terms  linear  in  ejit).  Since  superposition  does  not 
hold,  the  presence  of  e2  has  the  effect  of  "altering"  the  value  of  Klt  which  in 
turn  is  chosen  so  thatR(t)  is  minimized  in  the  RMS  sense.  In  general  the  out¬ 
put  y ( t)  of  the  nonlinear  element  consists  of  terms  involving  e^,  e2^,  and 
P  Q 

ej  e2  where  p  and  q  are  integers. 

The  observation  that  the  characteristic  N(e)  of  a  nonlinear  element  is, 
in  effect,  "altered"  by  addition  of  a  secondary  high  frequency  sinusoidal  signal 
e2(t)  to  the  primary  component  et(t)  was  observed  some  time  ago.  In  1945 
McColl  [1]  noted  the  apparent  linearizing  effect  of  the  "extra  oscillation"  on  a 
relay  control  system.  He  also  mentioned  the  idea  of  replacing  the  original 
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relay  element  with  a  new  "effectively  linear"  element  and,  by  so  doing,  avoid¬ 


ing  all  explicit  mention  of  the  extra  oscillation.  McColl  attributes  this  latter 
idea  to  Lozier  [2].  These  two  ideas  are  the  basic  notions  of  the  DIDF 
♦echnique;  they  actually  had  their  beginning  prior  to  any  published  work  on  the 
"single  input"  DF. 

Some  of  the  earliest  studies  of  describing  function  techniques  were 
conducted  independently  by  Goldfarb  [3],  Tustin  [4],  and  Kochenburger  [5]. 
These  early  researchers  were  evidently  inspired  by  the  so-called  method  of 
"harmonic  balance"  developed  somewhat  earlier  by  Kryloff  and  Bogoliuboff  [6]. 
The  success  and  wide  acceptance  of  both  the  DF  and  DIDF  can  be  attributed  to 
the  good  low  pass  filtering  characteristics  of  most  practical  controlled  plants. 
This  is  true  even  though  the  definition  of  both  the  DF  and  the  DIDF  depends  only 
upon  the  nonlinear  element  and  input  signals. 

The  use  of  a  triangle  wave  secondary  component  e2(t)  to  provide 
"precisely  linear"  signal  amplification,  in  the  sense  of  the  DIDF  for  carrier- 
controlled  relay  servos,  was  probably  first  proposed  by  Lozier  [2].  The  basic 
concepts  involved  in  the  DIDF  for  two  sinusoidal  components,  where  the  ratio 
of  frequencies  is  an  irrational  number,  were  clearly  outlined  by  Lozier.  West, 
Douce,  and  Lives ly  [7]  first  introduced  the  term  "Dual  Input  Describing  Func¬ 
tion"  and  made  a  major  contribution  toward  the  development  of  practical  com¬ 
putational  procedures.  However,  the  DIDF  of  West  et  al.  is  valid  only  for  two 
sinusoidal  components  where  the  frequency  ratio  of  the  two  sinusoids  is  a 
rational  number.  In  practice  this  frequency  ratio  is  usually  taken  as  an  integer 
or  the  reciprocal  of  an  integer.  The  DIDF  of  West  et  al.  is  very  cumbersome 
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because  of  the  many  parameters  involved.  In  a  later  study,  Oldenburger  [8] 
noted  the  stabilizing  effect  brought  about  by  the  addition  of  a  high  frequency 
sinusoidal  secondary  signal  into  the  nonlinear  element  input  in  certain  unstable 
nonlinear  feedback  control  systems.  Somnervil?'  and  Atherton  [9]  proposed 
a  very  appealing  approximate  method  for  calculating  multiple  input  describing 
functions  or  equivalent  gains  K  which  are  valid  for  several  input  components 
with  different  waveforms  including  dc,  sinusoidal,  and  gaussian  noise.  When 
the  input  consists  of  only  two  independent  components,  they  propose  and  outline 
a  method  for  obtaining  the  DIDF  by  a  two-stage  evaluation  process  in  which  the 
first  step  involves  the  calculation  of  an  effective  or  modified  nonlinear  element. 
Oldenburger  and  Boyer  [10]  utilized  a  modified  nonlinear  element  approach  to 
derive  many  useful  DIDF  curves  for  the  restricted  case  of  two  sinusoidal  com¬ 
ponents  where  the  frequency  of  the  secondary  component  is  high  compared  to 
the  primary  component.  The  modified  nonlinear  element  concept  was  later 
shown,  by  Gibson  and  Sridhar  [11 1,  to  give  exact  answers  (subject  of  course  to 
the  accuracy  of  computational  methods  and  equipment)  for  two  sine  wave  com¬ 
ponents  where  the  frequency  ratio  is  an  irrational  number.  The  DIDF  consid¬ 
ered  by  Gibson  and  Sridhar  was  derived  using  the  theory  of  random  functions 
and  is  valid  only  for  single  valued  nonlinear  elements.  Mu^h  earlier,  however, 
Bennett  [12] ,  and  later  Kalb  and  Bennett  [13] ,  used  a  double  Fourier  Series 
expansion  of  a  nonlinear  element  output  in  analyzing  modulation  products. 

Later  Amsler  and  Gorozdos  [14]  used  the  method  of  Bennett  and  Kalb  in  the 
analysis  of  bistable  control  systems  and  derived  what  was  later  called  the 
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DDIF  for  a  relay  with  two  sinusoidal  input  components  with  irrational  fre¬ 
quency  ratio. 

Other  investigators  such  as  Sridhar  and  Oldenburger  [15]  and 
Atherton  and  Turnbull  (16]  have  proposed  alternative  approaches  to  the  DIDF 
based  upon  statistical  methods.  Gelb  and  Vander  Velde  [17]  make  use  of  a 
simpler  dc  equivalent  gain  as  an  approximation  to  the  DIDF  for  two  sinusoidal 
components  in  the  analysis  of  limit  cycling  control  systems.  In  a  recent  book, 
Gelb  and  Vander  Velde  [18]  have  applied  the  integral  representation  of  non¬ 
linear  elements  used  by  Gibson  and  Sridhar  to  the  calculation  of  a  DIDF  for 
two  sinusoidal  input  components  (two-sinusoid  input  describing  function  or 
TSIDF) .  However,  in  their  book  Gelb  and  Vander  Velde  use  purely  determin¬ 
istic  arguments  rather  than  the  theory  of  random  functions.  In  this  same  book 
Gelb  and  Vander  Velde  also  disclose  a  very  appealing  power  series  TSIDF  cal¬ 
culation  for  odd  nonlinear  elements  with  input  sinusoids  whose  frequency  ratio 
is  irrational.  This  power  series  of  the  TSIDF  is  easily  generated  from  an 
algebraic  form  of  the  ordinary  DF. 

All  of  the  methods  discussed  above,  with  the  exception  of  the  DIDF  of 
West,  Douce,  and  Livesly  [7]  apply  only  to  single  valued  nonlinear  elements. 
Recently,  Mahalanabis  and  Nath  [19-21]  have  proposed  a  method  of  calculating 
the  DIDF  which,  they  assert,  holds  for  multivalued  nonlinear  elements  with 
multiple  inputs.  However,  Atherton  [22]  has  subsequently  pointed  out  that  this 
proposed  method  is  valid  only  when  the  nonlinear  element  may  be  described  by 
the  function  f(e,e) ,  where  e  is  the  input  to  the  nonlinear  element  [23].  Often 
a  function  f(e,e) ,  although  incorrectly  describing  the  characteristic  of  a  given 
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multivalued  nonlinear  element,  may  be  successfully  employed  when  the 
derivative  of  the  input  signal  does  not  change  signs  within  the  hysteresis  or 
multivalued  region.  This  condition  cannot,  in  general,  be  guaranteed  for  two 
independent  input  components  regardless  of  any  stipulations  on  the  amplitude  of 
these  signals.  Consequently,  the  method  proposed  by  Mahalanabis  and  Nath 
is  incorrect  for  two  input  components  even  when  the  restriction  is  added  that 
the  difference  in  amplitudes  of  the  two  components  must  be  greater  than  the 
hysteresis  width.  For  this  reason,  many  subsequent  articles  [24-28]  relating 
to  such  a  method  are  incorrect.  Xiohan  and  Krishma  [29]  have  correctly  used 
this  approach  to  find  DIDF’s  for  some  multivalued  nonlinear  elements  for  which 
the  function  f(e,e)  is  an  exact  representation.  This  rather  special  class  of 
nonlinear  elements  excludes  many  of  the  important  multivalued  nonlinear  ele¬ 
ments  encountered  in  practice,  however. 

Since  frequencies  and  the  frequency  ratio  of  the  two  periodic  components 
of  the  input  to  a  nonlinear  element  are  mentioned  quite  often  in  this  report, 
standard  symbols  will  be  adopted  for  these  terms.  Hereafter  the  frequency  of 
the  primary  input  component  (usually  a  sine  wave)  will  be  defined  as  co  and  the 
fundamental  frequency  of  the  deterministic  secondary  input  component  will  be 
defined  as  /?.  The  term  "frequency  ratio"  will  refer  to  the  ratio  =  /?/co. 

One  may  successfully  apply  the  modified  nonlinear  element  method  of 
determining  the  DIDF  for  multivalued  nonlinear  elements  when  certain  restric¬ 
tions  are  placed  on  the  input  components.  Atherton  and  Turnbuii  [16]  have 
demonstrated  that  good  approximate  answers  may  be  obtained  when  /?/co  »  1 
or  when/?/ co  «  1.  The  first  case  (/?/co  »  I)  requires  only  a  straightforward 
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application  of  the  method,  whereas  the  latter  case  ( {?/ «  1)  requires  a  slight 
modification  of  the  method  as  originally  outlined  by  Sommerville  and 
Atherton  [9]. 

Cook  [30]  has  given  the  DIDF  of  some  multivalued  nonlinear  elements 
for  both  sinusoidal  and  statistical  secondary  input  components  hy  using  the 
modified  nonlinear  element  method.  Cook  considered  only  the  more  straight¬ 
forward  application  of  the  modified  nonlinear  element  method  where  pj  >x  »  1. 

In  another  recent  paper  [31],  Mahalanabis  and  Oldenburger  have 
proposed  an  approximate  method  of  calculating  the  DIDF  of  a  multivalued  non¬ 
linear  element  by  the  use  of  statistical  methods.  They  assert  that  the  fre¬ 
quency  of  the  secondary  signal  component  may  be  either  higher  or  lower 
(p/u)  irrational)  than  that  of  the  regular  (primary)  input  component.  The 
latter  statement  conflicts  with  the  findings  of  this  study  and  will  be  discussed 
in  more  detail  in  Chapter  II. 

To  the  author's  knowledge ,  the  papers  cited  above  include  the  major  part  of 
the  published  original  work  on  the  DIDF's  of  multivalued  nonlinear  elements.  Only 
Mahalanabis  and  Oldenburger  [31  ]  assert  that  their  proposed  DIDF  applies  gener¬ 
ally  to  multivalued  nonlinear  elements  and  their  DIDF  appears  to  be  incorrect.  Even 
the  restrictive  case  where  the  nonlinear  element  input  consists  of  two  sinusoidal 
components  with  irrational  frequency  ratio  ($/ cehas  not  been  satisfactorily  solved. 

Seme  of  the  studies  cited  above  were  concerned  with  stability  consid¬ 
erations  and  others  were  concerned  with  the  signal  transmission  properties  of 
nonlinear  elements.  The  primary  concern  of  this  study  is  the  manner  in  whJch 
the  effective  nonlinear  characteristic  is  altered  in  the  presence  of  various 
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deterministic  secondary  signals.  The  possibility  of  changing  the  apparent  gain 
characteristics  of  a  nonlinear  element  by  injection  of  various  "stabilizing” 
signal  waveforms  has  been  considered  by  other  authors  [16,32,33].  Also  of 
concern  is  an  effective  method  of  dea:ing  with  multivalued  nonlinear  elements 
without  the  severe  restrictions  of  the  modified  nonlinear  element  concept. 

Since  it  is  known  that  the  injection  of  secondary  signals  of  different 
waveforms  at  the  input  of  a  nonlinear  element  results  in  different  DIDF's,  a 
related  synthesis  problem  may  be  posed.  This  synthesis  problem  is  stated  as 
follows:  Given  a  nonlinear  element  with  characteristic  N  (e) ,  find  a  waveform 
(if  one  exists)  of  a  periodic  secondary  input  component  e2(t)  which  will  result 
in  a  specified  DIDF.  With  added  qualifications,  this  problem  will  be  defined  as 
the  inverse  DIDF  problem.  Gibson,  Hill,  Ibrahim,  and  di  Tada  [34]  have 
proposed  an  inverse  DF  problem  where  it  is  desired  to  find  the  nonlinear 
characteristic  which  has  a  specified  describing  function.  Although  the  inverse 
DIDF  problem  is  not  a  logical  extension  of  the  inverse  DF  problem  defined  by 
Gibson  et  al. ,  it  is  perhaps  a  more  practical  one  for  the  tw:o-input  component 
case.  It  is  with  this  inverse  DIDF  problem  that  part  of  this  report  is  con¬ 
cerned.  An  unsuccessful  attempt  was  made  to  find  a  general  solution  to  this 
problem.  As  will  be  seen,  specific  classes  of  nonlinear  elements  lend  them¬ 
selves  to  relatively  simple  solutions.  In  the  absence  of  a  general  approach  to 

• 

the  inverse  DIDF  problem,  considerable  use  could  be  made  of  curves  showing 
the  DIDF's  of  the  nonlinear  element  for  several  specific  secondary  signal  wave¬ 
forms.  The  derivation  of  a  catalog  of  several  such  DIDF's  and  the  development 


of  new  and  shorter  methods  for  obtaining  them  comprises  another  contribution 
of  this  report.  The  specific  secondary  signals  considered  are  the  sine  wave, 
triangle  wave,  and  square  wave. 

In  summary,  it  is  felt  that  thus  report  makes  some  contribution  in  the 
following  specific  areas: 

a.  The  historical  aspects  of  the  DIDF 

b.  The  inverse  DIDF  problem 

c.  A  new  method  of  obtaining  the  DIDF  (The  proposed  method  holds 
for  a  broad  class  of  multivalued  nonlinear  elements  and  simplifies 
to  a  very  compact  form  for  single  valued  nonlinear 

elements  with  sinusoidal  input  comoonents. } 

d.  Calculation  of  DIDF's  for  several  specific  secondary  signals. 
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CHAPTER  II 

METHODS  OF  OBTAINING  THE  DIDF  FOR  DETERMINISTIC  INPUTS 


2.1  The  DIDF  of  West,  Douce,  and  Livesly 

The  most  general  DIDF  for  the  case  in  which  the  input  to  a  single  valued 
nonlinear  element  is  the  sum  of  two  sine  waves,  e  =  A  sin  (cat  +  <p) 

+  B  sin  (ncrt  +  ip) ,  is  given  by  the  expression 

1  2lr 

K(A,  B,n  ,  <p,4)  =  —  /  NfA  sin  {crt  +  <p)  +  B  sin  (neat  +  tp)  ]sin  (ut  +  4>)dtat. 

o 

(2.1) 

Such  a  DIDF1  has  live  variables,  A,  B,  n,  <p,  and  4>,  and  is  there  somewhat 
complicated  to  use  in  practical  problems.  West,  Douce,  and  Livesly  [7] 
developed  a  simplification  of  the  DIDF  given  by  Equation  (2. 1)  by  assuming  the 
parameter  n  to  be  a  rational  number.  In  fact  the  investigations  of  West  et  al. 
led  them  to  consider  the  even  more  restrictive  case  where  n  is  an  integer  or 


hereafter,  the  order  of  the  indicated  parameters  in  the  DIDF 
K(A,  3,a,b,n,  <p, . . . )  has  implied  meaning.  The  first  parameter  A  is  the 
amplitude  of  the  primary  input  signal  for  which  the  DIDF  is  derived.  The 
second,  B,  is  the  peak  value  of  the  secondary  input  signa1.  Then  follow  the 
parameters  a,  b,  c,  etc. ,  associated  with  the  nonlinearity  itself.  The  fourth 
group  (n,  4>,ip,etc. )  consists  of  the  frequency  ratio  and  phase  angles  associated 
with  the  two  input  signals.  These  comments  also  apply  to  the  modified  non¬ 
linearity  N(A,  B,a,b,c,n, . . , ) . 

13 


^£**2^^%****^™  -"7"'  — - ■“ - -7 — 


the  reciprocal  of  an  integer.  With  this  added  restriction,  the  inclusion  of  two 


phase  shift  quantities  <J>  and  ±  becomes  redundant  and  one  of  those  parameters 
was  eliminated.  In  this  way,  the  simplified  DIDF  of  W'est,  Douce,  and  Livesly 
was  obtained  in  the  form 

2jr 

*  1  c 

K(A,B,n,  (p)  =  — r  N[A  sin  (u>t  +  <p)  +  B  sin  not)  sin  (cut  +  $)dLjt  . 

7TA  j 

°  (2.2) 

when  N(e)  is  single  valued.  The  four  parameters  A,  B,  n,  and  <j>  in  the  DIDF 
of  West,  Douce,  and  Livesly  still  require  a  very  large  amount  of  data  to  give  a 
complete  representation.  For  this  reason,  and  because  of  the  restrictions  on 
n,  its  use  has  been  limited  primarily  to  investigations  of  stability,  subharmonic 
(superharmonic)  oscillations,  and  jump  phenomena  in  nonlinear  systems.  As 
pointed  out  by  West  et  al. ,  when  N(e)  is  adequately  described  by  a  low'  order 
polynomial  the  DIDF  is  most  easily  found  by  a  direct  expansion  to  obtain  the 
terms  in  the  output  with  frequency  ut.  Of  course,  the  direct  expansion  technique 
works  equally  well  when  the  frequency  ratios  of  the  input  sinusoidal  components 
are  irrational. 

2.  2  The  Modified  Nonlinear  Element  Concept 

The  concept  of  the  equivalent  nonlinear  element  or  the  modified  non¬ 
linear  element  probably  originated  with  Nikiforuk  and  West  f 35 ] .  However, 
their  modified,  normalized,  input-output  characteristic  was  defined  only  for  a 
sinusoidal  input  when  noise  was  added  to  this  input.  Sommerville  and 
Atherton  [y]  extended  this  concept  to  give  a  more  meaningful  "effective 
nonlinear  element.  "  This  effective  nonlinear  element  came  about  as  a  result 


of  the  two-stage  evaluation  of  the  equivalent  gain  of  a  nonlinear  element  with 
respect  to  a  primary  signal  in  the  presence  of  several  other  deterministic  or 
random  input  components.  The  restriction  was  imposed  that  the  cros'  correla¬ 
tion  function  of  any  two  of  the  random  input  components  must  be  zero  and  the 
frequency  ratio  for  periodic  input  components  must  be  irrational.  The  process 
is  shown  diagramatically  in  Figure  2.1. 

INPUTS 
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Figure  2. 1.  Diagram  Showing  the  Error  Term  R(x,y,z) 
and  Equivalent  Representation  of  a  Three-Input 
Nonlinear  Element  Where  R(x,y,z)  Is  Neglected 


In  general,  the  spectrum  of  the  output  V^ix.y.z)  consists  of  the  sum 
of  all  the  frequencies  contained  in  x,  y,  and  z  plus  harmonics  of  these  fre¬ 
quencies  and  frequencies  resulting  from  various  cross-products  of  x,  y,  and  z. 


The  equivalent  gain  K  is  that  value  of  K  (possibly  complex)  which  gives  a 

X  X 

minimum  value  of  R(x,y,z)  in  the  RMS  sense  as  is  allowed  to  vary.  The 

*  A 

other  equivalent  gains  and  K ^  are  found  independently  in  the  same  manner. 
Of  course,  if  only  the  equivalent  gain  to  one  input  component  (primary  input 
component)  is  desired,  as  is  the  case  in  this  study,  it  is  not  necessary  to  cal¬ 
culate  the  equivalent  gains  for  the  secondary  components.  When  the  primary 


input  component  x  is  a  dc  or  sinusoidal  signal  it  turns  out  that  the  required  out¬ 


put  K  x  is  simply  the  corresponding  dc  or  fundamental  Fourier  periodic  compo- 

X 

nents  of  V^x.y,  z)  at  the  frequency  of  x.  Sommerville  and  Atherton  [9]  have 

A 

shown  that  precisely  the  same  equivalent  gain  Kx  results  when  a  two-stage 
method  of  evaluation  is  used.  The  intermediate  step  is  to  define  an  effective 


nonliner  element  (modified  nonlinear  element)  by  considering  a  dc  signal  Aq 
instead  of  the  primary  signal  x  together  with  the  other  input  components.  The 

A 

function  N  ^A^  relating  the  average  dc  output  as  a  function  of  Aq  is  defined  as 
the  new  effective  nonlinear  element.  If  the  primary  input  component  x  is  then 

-A  A 

applied  to  the  new  effective  nonlinear  element  N,  the  signal  K^x  will  appear  at 
the  output.  This  concept  of  the  modified  nonlinear  element  gives  the  engineer  a 


very  helpful  physical  insight  into  the  mechanism  of  "signal  stabilization"  or 
equivalent  linearization  via  high  frequency  signal  injection  [33].  When  the 

A 

number  of  inputs  is  reduced  to  two,  the  equivalent  gains  K  are  called  DIDF's. 


As  mentioned  earlier,  when  the  input  consists  of  a  sum  of  sinusoidal 


f 


components,  an  equivalent  gain  or  DIDF  may  be  defined  with  respect  to  each 
component  in  the  input.  If  there  are  only  two  components,  both  of  which  are 
sinusoidal,  then  the  resulting  equivalent  gains  have  been  called  [36]  TSIDF's. 
Although  the  equivalent  gains  of  Somme rvi lie  and  Atherton  were  defined  and 
formulated  to  include  stochastic  components  in  the  input  signals,  their  formu¬ 
lation,  as  cited  in  this  report,  will  be  restricted  to  inputs  in  which  the  secon¬ 
dary  components  are  deterministic  periodic  signals.  Gelb  and  Vander 
Velde  138]  have  discussed  such  a  formulation  for  the  r,SIDF.  Suppose  the  two 
sinusoidal  components  of  the  input  are  given  by  B  sin  ;?t  and  A  sin  a)t,  where 
B/u)  is  an  irrational  number.  Then  the  equivalent  gain  or  TSIDF  with  respect 
to  the  component  A  sin  a)t  becomes 

1 

K(A,  B)  =  2tt^A  ^  sin  wt  dwt  J  N(A  sin  cut  +  B  sin  /St)d/3t  .  (2.3) 
o  o 

a 

A  similar  expression  defines  K(B,A)  with  respect  to  the  component 
B  sin  /3t.  It  should  be  pointed  out  that  Equation  (2.3)  holds  only  for  odd, 
single  valued,  nonlinear  elements.  As  remarked  earlier,  the  derivation  of 

a 

K(A,B)  may  be  carried  out  in  a  two-stage  process  by  first  defining  a  modified 

nonlinear  element  N^A^.B^  for  the  case  in  which  the  input  consists  of  the  sum 

of  a  dc  signal  A  and  a  sinusoid  B  sin  Bt.  The  characteristic  function  of  the 
o 

modified  nonlinear  element  then  becomes 

A  j  2  ir 

N(Ao*B)  =  /  N^o  +  B  sin/sAdiSt  ,  (2.4) 

and  the  TSIDF  is  given  by 
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(2.5) 


a  1  ^  A 

K(A,B)  =  — —  I  N(A  sin  wt,  B)  sin  cot  dut 
ttA  j 
0 

This  same  two-stage  process  will  be  used  in  this  thesis  except  that  the 
secondary  signal  is  net  restricted  to  be  sinusoidal.  Oldenburger  and 
Boyer  [10]  have  outlined  a  semi-intuitive  procedure  for  calculating  the 
TSIDF  w'hich  closely  parallels  the  modified  nonlinear  element  method  proposed 
by  Sommerville  and  Atherton.  In  the  outline  of  their  method,  Oldenburger  and 
Boyer  emphasize  the  concepts  which  underlie  the  experimental  measurement 
of  the  modified  nonlinear  clement  and  the  TSIDF. 

The  double  Fourier  Series  expansion  method  given  by  Gelb  and 
Vander  Velde  [18]  and  credited  to  Bennett  [12]  was  mentioned  in  Chapter  I. 
The  TSIDF  expression  given  by  this  method  is  exactly  that  of  Equation  (2.3) 
which  results  from  the  modified  nonlinear  element  method  and  will  not  be 
discussed  further. 

The  modified  nonlinear  element  method  given  by  Equations  (2.3) , 

(2.4) ,  and  (2.  5)  will  be  derived  in  section  2.6  of  this  chapter  by  an  approach 
which  differs  from  the  random  noise  approach  used  by  Sommerville  and 
Atherton  [9]  or  the  intuitive  approach  used  by  Oldenburger  and  Boyer  [10]. 


2.3  Power  Series  Method  (TSIDF) 

The  power  series  method  of  obtaining  the  TSIDF  mentioned  in 
Chapter  I  is  valid  for  odd  nonlinear  elements  when  the  frequency  ratio  of  the 
input  components  is  irrational.  This  method  is  the  result  of  expanding 
K(A,B)  in  a  Taylor  series  about  B  =  0.  Gelb  and  Vander  Velde  [18]  have 
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shown  that,  provided  ihe  necessary  derivatives  exist,  the  TSIDF  may  be 


written  as 


K(A.B)  =  ). - ^ -  V  (A)  ;  A  >  B 

p=0  (2)  *"P(  pi ):  P 


where  V  (A)  is  found  by  the  recursive  equation 

d2V  (A)  ,  dV  (A) 

V  (A1  P  ^  3  P 

p»l  ■  dA,  A  dA 


with  the  first  term  being. 


V  (A)  =  K(A) 
o 


the  ordinary  DF.  Likewise, 


(2.8) 


K(A,B)  =  y  — — - -  W  (E)  ;  A  <  B  , 

zn  n 


P=o  2”pp: (P  +  i)  p 


where 


d2W  <  B>  .  dW  (B) 

w  ,  B)  = — E — +  -i  - 

P"1  B  <iB 


with  the  first  term  being  found  by  the  relation 


(2.10) 


Wo(B)  -  K(B>  ♦  f  Sip 


(2.11) 


The  function  K(B)  is  the  ordinary  DF,  when  B  sin  fit  is  the  input  sinusoid.  The 
power  series  method  is  easy  to  apply  and  is  very  useful  when  closed  form 
solutions  for  the  TSIDF  cannot  be  found. 


2. 4  Methods  of  Stochastic  Processes 

Many  authors  [9,11,15,16,20,22,24,27,28,29,31,37]  have  used  the 


mathematics  of  random  functions  to  formulate  and  derive  equivalent  gains  and 
DIDF's.  Some  of  these  studies  differ  in  the  method  of  approach  but  all  involve 
the  same  fundamental  ideas.  Gibson  and  Sridhar  [111  used  this  approach  to 
derive  an  exact  DIDF  for  single  valued  nonlinear  elements  with  an  input  con¬ 
sisting  of  two  sine  wave  components  (TSIDF).  Mahalarabis  and  Nath  [37) 
proposed  a  more  direct  approach  to  calculating  the  TSIDF  via  the  techniques 
of  random  process  theory. 

The  TSIDF  for  single  valued  nonlinear  elements  was  shown  by  Gibson 
and  Sridhar  [111  and  Gelb  and  Vander  Velde  [18)  to  be 

QC 

K(A ,  E)  y  (ju)J0(Bu) .Jj(Au)du  ,  (2.12) 

_oo 

where 

oo 

y ( ju)  =  /  y(x)  exp[-jux)dx  (2.13) 

is  the  Fourier  Integral  Transform  of  the  nonlinear  element  output  y(x)  with  x 
as  the  input.  The  functions  J0(Ba)  and  Jj(Au)  are  Bessel  functions  of  order 
zero  and  one  respectively. 

2 . 5  Methods  of  Determining  the  DIDF  for  Multivalued 
Nonlinear  Elements 

As  already  discussed  in  Chapter  I,  relatively  few  results  have  beer: 
published  on  effective  methods  of  obtaining  the  minimum  RMS  error  DIDF  for 
multivalued  nonlinear  elements.  One  possible  exception  is  the  application 
[16,30,33]  of  the  modified  nonlinear  element  concept,  which  theoretically 
applies  [16  1  to  the  multivalued  nonlinear  element  case  only  as  n  or  1/n 
20 


E 


approaches  infinity.  It  is  plausible  that  the  DIDF  obtained  by  this  approach  is 
accurate  enough  for  practical  purposes  when  1  «  n  «  00 .  However,  to  the 
author's  knowledge,  no  previously  published  work  has  dealt  with  the 
inaccuracies  involved  in  this  finite  case  (1  «  n  «  00 )  or  with  the  question 
of  possible  ranges  of  finite  n  for  which  the  method  is  accurate  enough  for 
practical  purposes.  In  respect  to  the  application  of  the  modified  nonlinear 
element  method  to  multivalued  nonlinear  elements,  Cook  [30]  has  suggested  a 
method  of  keeping  track  of  the  phase.  A  new  approach  to  the  question  of  find¬ 
ing  the  DIDF  for  multivalued  nonlinear  elements  is  offered  in  the  following 
section.  This  approach  may  be  useful  in  determining  the  exact  DIDF  for 
multivalued  nonlinear  elements  with  deterministic  inputs. 


2 . 6  A  New  Average  DIDF  Method  for  Deterministic  Inputs 

Consider  the  general  TSIDF  given  by  Equation  (2,1)  and  let  n  be 
restricted  to  be  an  irrational  number.  This  is  the  case  considered  [9,10,11, 
16,25,26,27,31,32,33,36,37]  where  the  random  input  approach  was  used  to 
derive  the  TSIDF  for  deterministic  inputs.  In  such  applications,  the  random 
input  approach  merely  provides  a  convenient  method  for  averaging  out  the 
relative  phase  angles  <p  and  This  is  easily  understood  when  one  considers 
the  very  practical  problem  of  trying  to  experimentally  determine  the  TSIDF. 
When  n  is  irrational ,  the  phase  angles  <p  and  <P  lose  their  physical  importance 
and  it  becomes  impractical  to  even  include  them  in  the  TSIDF.  Consider  the 
accurate  experimental  setup  to  measure  the  TSIDF  given  by  Equation  (2.1). 

It  is  obvious  that  if  the  two  signal  generators  which  generate  the  components 
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A  sin  (cot  +  0)  and  Bin^t  +  0)  are  allowed  to  run  continuously,  a  different 
value  of  K  will  be  observed  with  each  succeeding  2r  interval  of  the  primary 

input  component  A  sin  (cot  +  0) .  What  is  actually  desired  is  the  average  value 
of  K  as  the  observation  time  increases  without  bound.  Therefore,  this 
"average  TSIDF"  will  not  include  0  and  o'  as  parameters.  In  other  words,  the 
actual  TSIDF  which  is  desired  is  given  by 

lim  1  2P~ 

K(A,B,n)  =  im  -  f  \'(Asincot  -»  B  sin  ncot ) sin  cot  dcot  (2.14) 

p— -00  prA  j 
o 

where  p  is  an  integer. 

Equation  (2.14)  is  equivalent,  in  the  measurement  problem,  to 
evaluating 

1  2? 

K(A,  B,n,  d>,  4)  =  —  I  N  [A  sin  (cot  +  4)  f  B  sin  (ncot  +  0)  )sin  (cot  +  4)  dcot 

it  J 

°  (2.15) 

for  all  possible  combinations  of  "phase"  or  zero  crossings  of  the  two  input 

A 

components  and  averaging  these  values  of  K  ( A ,  B,  n,  0,  4' ) . 

Thus,  there  is  another  way  to  calculate  the  TSIDF.  Let  one  of  the 
sinusoidal  input  components  contain  a  relative  phase  angle  (f>,  then  average  the 
resulting  set  of  "instantaneous  DIDF's"  as  0  ranges  over  2ir.  The  above 
arguments  are  more  rigorously  developed  in  the  ensuing  discussion.  Define 
an  equivalent  gain  gain  K  and  K  for  each  input  component  A  sin  wt  and 

n  D 

B  sin  (ncot  +  <p)  respectively.  The  "linearized"  output  is  given  by  the 

expression  K„  A  sin  cot  +  B  sin  (ncot.  +  d>) .  However,  K  and  K  are 
A  B  A  B 

complex  functions  of  A,  B,  and  n,  in  general,  and  can  therefore  be  expressed 
in  the  form 
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f 


K  =  K  +  jK  . 
A  ar  ai 


K. 


[B=Kbr+jKbi  '  (2’16> 

The  RMS  value  of  the  difference  between  the  actual  output  and  the  linearized 
output  over  all  time  is 

/  j.  j  2ttp 

R  -  I  *m  -  f  (N[A  sin  cot  +  B  sin  (ncot  +  d>)  ] 

\P—°°  2p-  ^ 


K  A  sin  cot  -  K  .  A  cos  cot 
ar  ai 


K,  B  sin  (ncot  +  d>) 
br 


2a.  \  % 


-  K,  .  B  cos  (ncot  +  <i>)}2dcot 
bi 


(2.17) 


where  p  is  an  integer.  Since  n  is  restricted  to  be  an  irrational  number , 
Equation  (2. 17 )  may  be  written  in  the  following  double  integral  form, 

III  2*  lTl 

R  =  ( —  —  /  d<p  J  (N[A  sin  cot  +  B  sin  (ncot  +  j 
\  '  o  o 


K  A  sin  cot  —  K  .  A  cos  cot 
ar  ai 


B  sin  (ncot  +  <p) 


B  cos  (ncot  +  *p)  ]  "dcot^  72  . 


(2.13) 


The  alue  of  the  real  part  of  K  (minimum  RMS  TSIDF1  which  minimizes  R  is 

A 


characterized  by  the  condition 


9R 


OK 


=  0  . 


(2.19) 


ar 
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The  indicated  differentiation  may  be  taken  under  the  integral  signs  since  the 

integrand  is  single  valued,  continuous,  and  differentiable  with  respect  to  the 

variable  K  even  though  N  ( e)  may  be  discontinuous  and  multivalued.  This 
ar 

a 

value  of  K  ,  defined  as  K  ,  is  found  to  be 
ar  ar 

A  1  27T  27T 

K  (A,  B,  n)  =  o*  f  dcp  f  N[A  sin  cut  +  B  sin  (ncut  +  <t>)  ]  sin  cut  dcut 
ar  Z7r  A  J  J 

0  0  „  „ 

(2.5 

Likewise  the  imaginary  part  of  the  DIDF  is  found  to  be 


K  .(A, B,n) 
ai 


2  2tt  2tt 

f  d<p  f  N[A  sin  cut  +  B  sin  (ncut  +  <p)  ]cos  cut  dcut 


2ttA 


(2.21) 


The  independent  nature  of  the  two  input  components  (n-irrational)  is  empha¬ 
sized  by  changing  the  variables  cut  and  ncut  to  Oi  and  02,  respectively. 

If  K  (e)  is  now  restricted  to  be  single  valued,  the  imaginary  part  K  . 

ai 

of  the  TSIDF  goes  to  zero  and  02  does  not  appear  in  the  answer  for  X  . 
Therefore,  Equation  (2.20)  may  be  reduced  to  the  following  form 
A  1  2ir  2ir 

K.  (A,  B)  =  0  2.  f  d<P  f  N(A  sin  6\  +  B  sin  <p) sin  Q<  dhi  .  (2. 22) 
A  Z7T  A  J  J 

0  o 

The  role  of  4>  and  02  tnay  also  be  interchanged  to  give  a  similar  expression  for 


the  DiDF, 


1  27T  27T 

Ka(A,  B)  =  2^2A  J  d^2  f  N(A  sin  0,  +  B  sin  02)sin  d 0,  .  (2. 23) 
oo 

.Since  the  integrand  is  single  valued  with  at  most  a  finite  number  of  finite 
discontinuities,  the  order  of  integration  may  be  interchanged.  That  is, 
Equation  (2.23)  may  be  written  as 
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1  27T  2tt 

Ka(A,  B)  =  2 A  j  sin  °i  <i61  J  N(A  sin  0j  +  B  sin  02)d#2  •  (2-24; 

o  o 

Equation  (2.24)  is  equivalent  to  (2.3)  which  in  turn  is  a  result  of  the  modified 
nonlinear  element  method.  Equation  (2.24)  is  also  the  form  of  the  TSIDF 
given  by  Gelb  and  Vander  Velde  [36  j .  The  two-step  evaluation  of  Equation 
(2.24)  may  be  carried  out  as  indicated  in  Equations  (2.4)  and  (2.5) . 

Although  the  order  of  integration  is  not  important,  the  integral  forms  resulting 
from  Equation  (2.23)  may  be  easier  to  integrate  than. those  of  Equation  (2.24) , 
or  vice  versa. 

Although  the  minimum  RMS  TSIDF  was  derived  above  under  the 
assumption  that  the  secondary  component  in  the  input  was  a  sine  wave, 

B  sin  (na,t  +  </>) ,  this  restriction  is  not  necessary.  In  fact,  any  periodic 
secondary  signal  cr{(3t)  could  have  been  chosen  without  affecting  the  form  of 
the  answer.  This  result  will  be  used  in  Chapter  V  to  show  how  the  DIDF  is 
affected  by  changing  the  wave  shape  of  the  periodic  signal  cr  (/3t) .  Equations 
(2.22)  through  (2.24)  also  emphasize  that  the  TSIDF  of  single  valued  non¬ 
linear  elements  is  independent  of  the  frequency  ratio  n.  It  appears  that  such 
a  simplification  cannot  be  made  when  considering  multivalued  nonlinear 
elements.  This  latter  statement  is  not  in  agreement  with  the  assertion 
recently  made  by  Mahalanabis  and  Oldenburger  [31] .  Further  discussion  of 
this  point  is  made  in  Chapter  IV  where  the  DIDF  for  the  relay  with  hysteresis 
is  found  by  use  of  Equations  (2.20)  and  (2.21)  and  is  compared  with  the 
published  results  of  Mahalanabis  and  Oldenburger. 
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Since  the  frequency  ratio  n  dees  not  enter  into  the  TSIDF  for  single 
valued  nonlinear  elements,  Equation  (2.20)  suggests  an  alternative  approach  to 
the  calculation  of  the  TSIDF.  For  instance,  when  n  is  very  close  to  one,  the 
DIDF  may  be  written  in  the  form 


2tt  2tt 


K(A,B)  =  ^  -  f  d <p  J  N[Dsin(cut  +  <£)  ]sin  cut  dcut  ,  (2.25) 


where 


D  =  D(A,  B,  4))  =  n/  A2  +  B2  +  2AB  cos  0 

/  .  _i  /  B  sin  <b 

sm  7  ? 


(2.26) 


;  B  cos  <p  +  A  2-  0 


7 r  -  sin 


\^J  A''  +  B  +  2AB  cos  (p 
B  sin  <f> 


(2.27) 


ln/a2 


1;  B  cos  <p  +  A  £  0  . 


+  B  +  2AB  cos  cf)J 
Now  evaluate  K(A ,  B)  by  a  two-stage  process,  similar  to  that  used  in  the 
modified  nonlinear  element  concept,  where 


g(A 


■ B •*>  ‘ikf 


N[D  sin  (cut  +  9)  ]sin  cut  dcut 


(2.28) 


and 


1  2n 

K(A.B)  =—  f  g(A,B,4>)d4>  .  (2.29) 

27T  J 
0 

The  function  g(A,  B,  (p)  simply  defines  a  special  case  of  the  DIDF  of 

A 

West,  Douce,  and  Livesly  where  n  =  1.  Equation  (2.29)  shows  that  K(A,B)  is 
the  average  of  g(A,B,<p)  as  <p  ranges  over  27r.  This  implies  that  the  average 
value  of  the  DIDF  of  West,  Douce,  and  Livesly  as  $  ranges  over  (0,  2tt)  is  the 
same  regardless  of  the  value  of  n.  The  function  g(A,  B,  <j>)  has  the  appearance 
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of  the  conventional  DF  with  the  additional  parameter  3>.  It  will  be  shown  oelow 
that  the  computation  oi'  g(A,B,  <*>)  requires  the  same  effort  as  the  computation 
of  the  DF.  It  may  be  considerably  more  difficult  to  perform  the  integration 
indicated  in  Equation  (2.29)  than  that  in  Equation  (2.  28) .  Since  the  charac¬ 
teristic  I\<e)  is  single  valued,  the  phase  0  is  preserved  through  the  nonlinear 
operation  indicated  in  Equation  (2, 28) .  Therefore,  Equation  (2.28)  maybe 
rewritten  as 


,  2- 
cos  v  r 

g(A,B,0)  = - - —  N(D  sin  cut)  sin  cut  dcut 

7rA  j 


(2.30) 


Substituting  Equation  (2.27)  in  Equation  (2.30)  gives 

27 r 


g(A,  B,  <p)  -- 


A  +  B  cos  0> 


3_ 

n/  A2  +  B2  +  2AB  cos  4 >  o 


J  N  (D  sin  cut)  sin  cut  dcut 


(2.31) 


and 


K(A,  B) 


277 


(1  +  B/A  ccs  4>)dip 
A2  +  B2  +  2AB  cos  (b 


2-  _ 

f  N (\!  A2  +  B2  +  2AB  cos  <p  •  sin  cut)  sin  cut  dcut  .  (2.32) 

o 

It  is  interesting  to  note  that  this  formulation  leads  to  a  general  expression, 
Equation  (2.32)  ,  which  has  the  appearance  of  an  elliptic  integral.  However,  it 
appears  that  this  integral  may  be  difficult  or  impossible  to  reduce  to  one  of  the 
standard  or  tabulated  elliptic  integrals.  The  functional  cos  [0(A,B,0)  ]  does 
not  depend  upon  the  nonlinear  element  and  is  symmetrical  about  <p  =  tt.  In  fact 
g(A,B ,<p)  is  symmetrical  about  <p  =  tt  as  well.  In  some  cases  approximate 
algebraic  solutions  may  be  desired  rather  than  graphs  of  the  TSIDF.  If 
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trapezoidal  integration  or  other  approximate  techniques  are  used,  accuracy 
sufficient  for  practical  application  is  usually  achievable  with  very  few  terms. 

0>B)  is 

oridinarly  not  difficult  to  obtain.  The  main  task  consists  of  determining  the 
fundamental  component  of  the  output  when  the  primary  component  of  the  input 
is  operated  upon  by  the  function  N  ^Aq,  . 

Example:  DIDF  of  Relay 

In  order  to  demonstrate  the  utility  of  the  average  DIDF  method,  con¬ 
sider  the  DIDF  of  the  perfect  relay.  The  relay  is  described  by  the  equation: 


In  the  modified  nonlinear  element  method,  the  function  N  (a 


y  =  M  sgn(e) 


M;  e  >  0 
-M;  e  <  0  , 


(2.33) 


where,  for  the  present  case,  e  =  A  sin  cot  +  B  sin  /3 t.  Equation  (2.30)  is 
simply  the  conventional  DF  times  cos  4>  with  D  sin  cot  as  the  input.  Thus, 


27 r 


1  r 

DF  =  — —  I  M  sgn(D  sin  cot)  sin  cot  dcot  = 
7rA  j 


4M 

7tA 


(2.34) 


and 


.  .  _  . ,  4M  /  1  +  B/A  cos  <p 

g( A,  B,  <p)  =  —  (  - - - - 


11  Vi  A2  +  B2  +  2AB  cos  <p  ) 
The  second  integration  indicated  by  Equation  (2.29)  gives 

A  4M  f  1  +  B/A  cos  i>  , , 

K(A,B)  =— -J  _  d <p 

o  n/  A2  +  B2  +  2AB  cos  0 


(2.35) 


(2.36) 


HtSHW'' 


f 


This  integral  may  be  evaluated  [38]. 


K(A,B)  = 


7T2  A  + 


H[(1*T£P1)p‘i>**TFE<w]  • 


(2.37) 


where  F(k)  is  the  complete  elliptic  integral  of  the  first  kind  and  E(k)  is  the 


complete  elliptic  integral  of  the  second  kind.  The  modulus  k  is  given  by 

;  b/A 


1  +  B/A  ‘ 


(2.38) 


Equation  (2.37)  is  valid  for  all  A  and  B.  However,  a  modulus  transformation 
may  be  applied  [39]  whereby 


2E(k)  -  (1  -  k2)  F(k) 


1  ^  k 


(2.39) 


F(k)  =  (1  +  k)F(k)  , 


(2.40) 


when 


k  -  B/A  ;  B  <  A 


k  =  A/B  ;  B  >  A 


(2.41) 


When  Equations  (2.38)  through  (2.41)  are  substituted  in  Equation  (2.37) ,  the 
following  expressions  for  the  DIDF  of  the  perfect  relay  results: 


K(A,B) 


8m 

—  E(B/A)  ;  B  <  A 

7T 


8M  b  F  /  a2\ 

^-|E(A/B)  -  (l  -  ^-jK(A/B)J  ;  A  <  B  .(2.42) 


Equation  (2.42)  coincides  exactly  with  the  expressions  obtained  by  Gibson  and 
Sridhar  [11]  and  Mahalanabis  and  Nath  [37],  who  used  essentially  different 


methods  of  calculations. 

vVhen  the  secondary  input  component  is  nonsinusoidal,  a  simplified  form  of 


the  average  DIDF  method  comparable  to  that  given  in  Equation  (2.32)  cannot  be 


found.  However,  either  of  the  expressions  given  in  Equations  (2.23)  and 
(2.24)  may  be  used  when  N(e)  is  single  valued.  Expressions  similar  to 
Equations  (2.20)  and  (2.21)  are  used  whenN'(e)  is  a  multivalued  function.  In 
order  to  include  a  more  general  periodic  secondary  signal  cr(/?t) ,  Equation 
(2.23)  can  be  written  as  follows: 


2  2ir  ?~ 

K[A:c\-,)]  =  ~~J~  f  d02  f X[A  sin  0  -  o( 02)  ]sin  0,  dft,  .(2.43) 
2tt  A 

0  o 

The  evaluation  of  Equation  (2.43)  may  be  carried  out  in  two  steps,  following 
the  modified  nonlinear  element  approach,  to  give 


‘if/1 s[a„  * 


do. 


(2.44) 


K[A;  a(  nz)  ) 


N[A  sin  (1,  ; 


cr  (02)  ]dfl, 


(2.45) 
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CHAPTER  ID 

DEFINITION  AND  SOLUTION  OF  AN  INVERSE  DIDF  PROBLEM 
FOR  A  PARTICULAR  CLASS  OF  NONLINEAR  ELEMENTS 


Although  a  large  portion  of  this  report  is  devoted  to  the  analysis  of  the 
output  of  a  nonlinear  element  when  the  input  contains  a  secondary  signal  of 
various  waveforms,  the  motivation  for  this  study  was  provided  by  the  somewhat 
more  difficult  synthesis  problem  of  finding  the  secondary  signal  o(/3t)  which  will 

A 

produce  a  prescribed  modified  nonlinear  function  N(A,  B)  for  a  given  nonlinear 
function  N(e) .  It  is  assumed  that  the  given  nonlinear  element  is  an  inherent 
part  of  a  system  and  cannot  itself  be  changed  or  replaced.  Perhaps  the  first 
logical  step  in  a  synthesis  procedure  is  to  define  the  desired  DIDF 

A  A 

K  =  K[A;  v([3t)  ]  where  the  primary  component  of  the  input  signal  is  a  sine 
wave.  If  one  then  attempts  to  work  backward,  using  the  modified  nonlinear 
element  concept,  the  second  step  consists  of  finding  the  equivalent  nonlinear 

A  A 

element  N  =  NjA  ;  cr(/?t)J  which  gives  the  desired  DIDF.  This,  or  course, 
is  exactly  analogous  to  the  inverse  describing  function  problem  posed  and 
largely  solved  by  Gibson,  Hill,  Ibrahim,  and  di  Tada  [34].  Even  though  a 
method  exists  which  may  be  used  to  solve  this  second  step,  one  might  proceed 
alternatively  by  first  defining  the  desired  nTa^;  <r(j3t)l  rather  than  the  desired 


t 

k 

% 


i 


K[A;  cm/? t)  j.  This  latter  problem  is  important  in  its  own  right  and,  in  most 
practical  situations,  might  just  as  well  be  the  starting  point.  A  logical 
extension  of  the  ideas  of  the  inverse  DF  problem  mentioned  above  would  lead 
one  to  define  the  inverse  DIDF  problem  as  one  of  finding  a  nonlinear  element 
which  corresponds  to  a  given  TSIDF.  However,  it  is  felt  that  the  inverse 
DIDF  definition  given  here  is  probably  a  more  practical  one.  The  inverse 
DIDF  problem  to  be  considered  here  may  be  precisely  stated  as  follows: 

The  Inverse  DIDF  Problem 

Given:  A  time  invariant  nonlinear  element  with  characteristic 
N(e)  and  the  desired  modified  nonlinear  characteristic 
N*  =  nJAq;  or(/3t)J 

Find:  The  periodic  secondary  signal  <7*  (/5t)  which  will  produce 
the  modified  nonlinear  characteristic  N*  . 

The  diagram  of  Figure  3. 1  shows  the  essential  elements  and  signals  required 
in  the  development. 

In  the  present  discussion,  only  odd,  single  valued  nonlinear  elements 
whose  characteristic  is  not  a  function  of  frequency  will  be  considered.  In 
order  to  conveniently  choose  and  restrict  the  class  of  secondary  signals  which 
should  be  considered,  the  following  observation  on  integrals  of  functions  is 
made:  Let  f(t)  be  a  given  function  of  time,  Aq  be  an  arbitrary  constant,  and 
N(e)  be  a  time  invariant  single  valued  nonlinear  function.  In  addition,  let 
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(a)  Exact  Representation 


DESIRED  MODIFIED 
NONLINEAR  ELEMENT 


(b)  Equivalent  Representation 

Figure  3.1.  Diagram  Illustrating  the  Inverse  DIDF  Problem 


Then  there  exists  a  function  g(t)  such  that 


dg(t)  >  0  :  0  <  t  <  t  , 

dt 

(3.2) 

and 

ti 

r  v 

J  N|A  +  g(t)  dt  =  X  . 

(3.3) 

o  L  °  J 

Therefore,  following  the  definition  of  the  modified  nonlinear  element  given  by 
Equation  (2.44) ,  one  may  consider  only  secondary  components  a(/3t)  belonging 
to  the  class  Z  where 


I 
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cr(,3t)  =  -a  (-fit) 


Yj :  cr(*/2  +  t)  -  ct(7t/2  -  /3t) 


da(>lt) 


>  0  ;  0  £  pt  *'  tt/2 


This  restriction  does  not  eliminate  from  consideration  any  modified  charac¬ 
teristic  function  N  |a  ;  cr(/3t)J  which  might  otherwise  be  considered. 

The  restriction  to  the  class  Z  greatly  simplifies  the  search  for  a  solu¬ 
tion  ct*  (/3t)  and  does  not  eliminate  any  useful  solutions.  It  is  obvious  that  a 
secondary  signal  can  alter  the  modified  nonlinear  characteristic  only  over 
certain  regions.  Such  aregion  in  the  Nf"A  ;cr(/?t)J  versus  plane  might  be 

as  shown  in  Figure  3.2.  In  addition  there  will  be  limits  on  the  derivatives  of 

,m  » 

A  d  X 

N  in  this  region  R.  That  is.  the  values  of  the  derivatives  - -  cannot  be 


arbitrarily  specified.  To  define  the  region  R  and  the  limits  on  the  derivatives 


which  may  be  achieved  by  a  secondary  component  in  the  input  is  a  large 


task.  However,  with  the  exception  of  specifying  the  region  R  a  priori,  this  is 
what  is  required  in  the  inverse  DIDF  problem.  Curves  showing  the  DIDF's  of 
nonlinear  elements  for  several  specific  secondary  signal  waveforms  should 
lend  some  insight  into  what  may  be  expected  from  the  inverse  problem.  A 
catalog  of  several  such  DIDF's  is  given  in  Chapter  IV. 

Consider  the  class  of  nonlinear  elements  defined  by  the  equation 


N(e)  =  Ke  sgn  (e)  (3. 

where  m  is  a  non-negative  even  integer.  Let  the  class  of  desired  equivalent 
nonlinear  elements,  N*  ,  be  defined  by  the  relation 


w 


I 


ft  [*,;.(«] 


Figure  3.2.  Region  over  Which  N  |a  ;  <r(/?t)J 
May  Be  Varied  by  Changing  cr  (/3t) 

N*  =  IsTA  ;<r*(/?t)l  -  k,  A  n  sgn  (a  \  .  (3.6) 

L  o  J  o  \  o/ 

This  equation  holds  for  some  as  yet  undefined  region  about  the  origin  A^  =  0. 
The  equation  for  the  equivalent  nonlinear  element  with  any  secondary  signal  and 
a  single  valued  nonlinear  element  is  given  by  the  functional 

1  2” 

nJaq;  a(/?t)]  =  —  /  +  o  (j3t)J  d/?t  .  (3.7) 


The  class  of  problems  under  consideration  can  be  solved  with  some  degree  of 

th 

success  in  a  very  straightforward  manner  by  defining  only  the  m  derivative 
of  N*  .  Of  the  nonlinear  elements  described  by  Equation  (3.5),  only  the 
perfect  relay  and  absquare  will  be  considered.  However,  the  method  proposed 
here  may  be  carried  further  witnout  defining  the  nonlinear  element  to  be  con¬ 
sidered.  Substitution  of  Equation  (3.5)  into  (3.7)  results  in  the  following 
desired  modified  nonlinear  element: 


[Ao;  cr*  (/?t)j  =  —  /  [Aq  +  a*(0t)jmsgn|^Ao  +  <r*  (£t)J  d/5t  .  (3.8) 
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It  is  necessary  to  separate  the  integral  of  Equation  (3.8)  into  two 


e 
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It  is  evident  from  Equations  (3.9)  and  (3.14)  that 


d  N:< 


>  0  . 


(3.18) 


The  latter  relation,  in  conjunction  with  Equation  (3. 11)  ,  requires  that 


n  i  m  -  1 


(3.19) 


except  when  m  =  0.  Therefore,  the  order  of  the  equation  describing  the  modi¬ 
fied  nonlinear  element  cannot  be  made  lever  than  m  -  1.  However,  the  order 
of  the  equation  describing  n|a  ;  cr </3t) J  may  be  larger  than  m.  Equation  (3. 10) 
describes  the  average  output,  of  an  absolute  value  circuit  and  is  always  positive 
when  and  a(/3t)  are  not  identically  zero.  Thus,  the  constant  of  Equation 
(3. 12)  is  nonzero  except  possibly  when  n  =  m  -  1.  The  proposed  method  is 
further  illustrated  by  the  following  examples: 

Example  1  (Perfect  Relay) 

til  A 

In  this  case,  m  =  0  and  only  the  0  derivative  of  N  is  required  to  find 
the  secondary  signal.  From  Equation  (3.17)  the  required  secondary  signal  is 


1  1 


-»-(t ;)'(£)".  •' 


j?t  <  tt/2  . 


(3.20) 


When  n  =  1,  that  is  when  the  relay  is  "linearized,"  the  signal  <j*  (/St)  is  a 
triangle  wave.  The  range  over  which  Equation  (3.2)  holds  is  limited  by  the 
amplitude  B*  of  the  secondary  signal.  The  amplitude  B*  for  the  relay  case  is 
readily  seen  from  Equation  (3.  20)  to  be 


/k_Y 

[kj  ’ 


(3.21) 
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The  two  degenerate  cases  where  n  =  0  anri  l-  -  n 

wucre  n  u  and  k,  -  0  must  be  considered 

separately.  When  n  =  o,  N*  =  k.  <?<m  / a  \  u-  u  ^ 

Sgn  (Ao)  ’  whlch  describes  the  original 

non, inear  element  and  .  «  „  the  only  correct  solution  „  „  js  ^ 

'■o  make  k,  -  0.  Epuation  ,eads  t0  a  wave  wlt„  M„ite  ^ 

by  referring  to  Equation  (3.15)  one  obtains  the  relation 

2K01  _  . 

*  kl  =  0  ’  Ac  <  B*  •  (3.22) 

The  last  equation  implies  that  =  0  for  nonzero  values  of  Aq.  One  signal 

0  W  SallSfylng  this  eQuation  is  a  symmetric  square  wave  of  arbitrary 
amplitude  B* .  The  resulting  N  is,  of  course,  a  relay  with  dead  band  width 

2B*  [161  *  N°tiCe  th"  3  "  “  *  negative  §1.  Solving  Equation 

(3.8)  when  m  =  0  gJves  ° 


K  ^7r 

N*  =  2^J  Sgn  [A0  +  o*(/?t)l  d/?t  =JHgt  . 

O  L  J  7T 


Thus 


(3.23) 


dN*  =  2K  d 0j_ 

dA  7r  dA  ’ 
0  o 


which  cannot  be 


(3.24) 


“CSa  ve.  as  is  easily  seen  in  Figure  3.3.  The  inequality 
n  *  0  must  therefore  be  Imposed  in  the  relay  problem.  A  typical  plot  of  N» 
for  n  =  1/2  and  k,  =  K  is  shown  in  Figure  3.4. 

Example  2  (Absquare) 

The  solution  for  the  abspuare  problem  ,s  obtained  by  simply  setting 

m  =  2  EqUaH0"  W 1?) '  The  ~  **  *■»  1/4  wave  of  the  secondary 

signal  is 


4 

M 
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Figure  3.4.  Modified  Nonlinear  Element  and  Secondary  Signal 
for  Perfect  Relay  Where  N*  is  Given  by  kj 
over  a  Positive  Range  of  A 


am  = 


2K 
ktn(n  - 


1) 


n-2  /  ,n-2 

p  l 


-777  ;  0  ^  pt  <  tt/2  .  (3.25) 


By  specifying  the  second  derivative  to  be 


d2N*  ,  , , ,  ,  A  n-2  , 

dX"T  =  n(n  -  l)kj  iAq  lsgnAQ 


(3.26) 


it  is  only  possible  to  guarantee  that 


dN*  ,  ,  n-1 .  _ 

~77—  =  nkJA  Isgn  A  +  C, 
dA  1  o  o  1 

o 


(3.27) 
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and 

X*  -  kjlA^lsgn  ^A^  +  C,Ao  -  C2  ,  (3.28) 

where  C,  and  C2  are  constants.  However,  solution  of  Equation  (J.  8)  when 
A^  =  0  shows  that  N[0,  <r*(f>t)]  =  0  for  allowable  a*  (fit),  and  thus  C2  =  0. 

Three  degenerate  cases  occur  when  n  =  0,  1,  2.  The  first  one,  n  =  0, 
is  eliminated  by  Equation  (3. 19) .  When  n  =  1,  that  is  when  the  absquare  is 
linearized,  it  is  convenient  to  use  Equation  (3. 15)  from  which  one  may  write 

0,  =  0  ;  Aq  *  °  .  (3.29) 

In  addition,  Equation  (3.13)  requires  that 

o(0)  =  Aq  ;  Aq  <  M  ,  (3.30) 

where  Mis  arbitrary.  The  only  secondary  signal  which  satisfies  the  require¬ 
ments  of  Equations  (3.29)  and  (3.30)  is  one  with  infinite  slope  for  cr (/3t)  <  M. 
For  example,  the  secondary  signal  of  Figure  3.5  results  in  the  modified  non¬ 
linear  element  of  Equation  (3.31)  which  is  linear  in  A  . 

o 


Figure  3.  5.  A  Secondary  Signal  Which  Linearizes  the  Absquare 
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n[a  ;  (T*  (/?t)l  -  Kr2M  +  i=ilA  ;  A 

L  O  J  L  7TJ0  O 


<  M  . 


(3.31) 


The  family  of  curves  for  N  and  K  are  given  in  Chapter  IV  when  Bj  =  0  and  M 

takes  on  different  values.  As  stated  above,  the  constant  Cj  of  Equation  (3.28) 

is  not  zero  for  values  of  n  other  than  one.  The  case  where  n  =  2  describes 

the  original  nonlinear  element  and  cr (/3t)  =  0  is  the  only  correct  solution.  The 
1 

exponent  — - -in  Equation  (3.5)  causes  the  secondary  signal  to  have  infinite 

amplitude  when  n  <  2.  Therefore  the  specified  N  cannot  be  physically 
realized  for  the  admissible  ramge  1  <  n  <  2. 

The  last  case  which  will  be  considered  under  the  absquare  example 
occurs  when  n  =  3.  The  secondary  signal  becomes 


•«*»  ■  0  ***«'*  ■ 


(3.32) 


which  is  a  triangle  wave  with  amplitude  B*  =  — — .  Substituting  this  wave  for 

i)  Kj 

a 

a(j3t)  in  Equation  ,(3.23)  jives  the  following  equations  for  N: 

n[V’H  -  K[siAoJ  *  8,1  Ao]  ;  \  “  84  (3-33> 

n|Aq;  <x(j3t)J  =  K^AQ2  +  yB*2J  ;  A0  -  B*  •  (3.34) 

K 

Therefore  the  gain  factor  kj  in  Equation  (3.27)  becomes and  the  constant 
Ct  becomes  KB* .  Although  the  kjA^3  term  does  appear  in  the  expression  for 
N  A  ;  a(j3t)  ,  it  is  not  large  compared  with  the  CjA  term  except  when  B  «  1, 

L  o  J  0 

in  which  case  the  usable  range  is  much  less  than  one.  The  results  in  this  case 
are  not  very  satisfactory.  It  is  suspected,  however,  that  not  any  other 

A 

secondary  signal  will  give  an  N  closer  to  the  one  specified. 


42 


» 


Consider  next  the  nonlinear  elements  which  are  defined  by 


N(e)  =  KeF  , 


(3.35) 


where  p  is  an  odd  positive  integer .  It  will  be  shown  that  all  secondary  signals 

a 

give  the  same  structural  form  of  N.  Differences  in  secondary  signals  alter 

/K 

only  the  magnitude  of  certain  parameters  in  N.  The  modified  nonlinear 
element  is  given  by 


1 

N[Ao;  <r(/3t)]  =  —  /  K [Aq  +  a(/3t)] Pc^t  .  (3.36) 


The  p  and  lower  order  derivatives  are  computed  as  follows: 


d'N  1  r 

— -r p:  2 r/Kd(5t  =  p:K 

dA  o 

o 


(3.37) 


d^N 


j£T  *  p'-  IT-f  [Ao  +  p<w]  dM  "  plKA  (3-33) 


o 


where 


[V+C>]  • 


1  Z7r 

c'  =  27  /[CT<w]  V 


=P!K/A3  +  3AC\ 

dA  P'3  6  '  0  0  ' 

o 


and  so  on  for  the  lower  order  derivatives. 


(3.39) 


(3.40) 


(3.41) 
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When  <r(/3t)  ^  0,  all  the  constants  which  appear  are  of  the  form 


1  9 

cn  =  27/  Ffr(^t)  ]  nd.^t  (3.42) 

0 

and  are  nonzero.  For  example,  the  modified  nonlinear  element  for  the  cubic 
is 

n[a  ;<r(/Jt)]  =  Aq3  +  3A  .  (3.43) 

Thus  by  changing  the  waveform  of  cr(j3t)  only  the  magnitude  of  Cj  is 
affected.  It  is  found  in  practice  that  the  equivalent  dc  gain  K  [a^;  cr(|3t)J  is 
often  quite  similar  to  the  DIDF  and  can  be  obtained  by  dividing  n|a  ;a(/?t)J  by 
Aq.  The  dc  equivalent  gain  for  the  cubic  is  given  by  the  relation 

K[Ao;  cr(/3t)J  -  Aq2  +  3Cj  .  (3.44) 

Thus  the  gain  for  small  inputs  ^Aq  —  0^  is  changed  from  zero  to  3Cj.  This  is 
a  peculiar  class  of  nonlinear  elements  where  the  waveform  of  cr  (/3t)  does  not 
affect  the  "shape''  of  N.  Therefore  the  secondary  signal  is  not  unique  even 
within  the  restricted  class  of  signals  described  by  Equation  (3.4) .  The  values 
for  the  constant  Cj  when  u(/?t)  is  a  square  wave,  sine  wave,  and  triangle  wave, 
respectively,  are 


(square  wave)  =  B2 

(3.45) 

Cj  (sine  wave)  --  B2/2 

(3.46) 

Cj  (triangle  wave)  =  B2/3  . 

(3.47) 

The  method  of  taking  derivatives  of  N  to  get  simpler  expressions  with 
which  to  work  is  not  limited  entirely  to  the  nonlinear  elements  defined  by 
Equation  (3.5) .  Limited  success  may  be  achieved  where  the  nonlinear 
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element  is  piecewise  constant  or  linear.  For  example,  the  "linearization"  of 
the  relay  with  dead  band  by  the  method  presented  above  requires  a  triangle 
wave.  The  triangle  wave  is  known  to  be  a  correct  answer  'Chapter  IV) .  The 
method  also  gives  the  same  answer  for  the  relay  with  dead  band  and  the 
perfect  relay  when  n  =  2.  The  secondary  signal  is  computed  to  be 

V  V 

<tW)  =  {y-)  (^2)  ;  0  ^  0t  <  tt/2  .  (3.48) 

However,  the  process  is  much  more  tedious  in  the  case  of  the  relay 

* 

with  dead  band  and  the  resulting  N  is  not  what  one  might  expect.  The  specified 
structure  |^N  =  kjA^2  sgn  ^A_^j  is  achieved  (within  a  constant)  only  in  the 
region  a  ^  !Aq!  s  B  -  a.  A  linear  range  results  when  I Aq  I  <  a 
(Figure  3.6) . 

An  important  question  which  has  been  left  unanswered  by  this  procedure 
is:  Does  a  secondary  signal  exist  which  gives  the  desired  modified  nonlinear 
element  J^N  =  Aj  sgn  ^A^J  ?  The  procedure  used  above  also  suffers  from  the 
fact  that  it  cannot  be  used  when  a  more  general  type  nonlinear  element  is 
encountered. 

A  brute  force  method  which  offers  some  possibility  for  the  general 
problem  is  now  given.  Only  odd  single  valued  nonlinear  elements  will  be  con¬ 
sidered.  The  method  is  approximate  and  involves  assuming  the  amplitude  B 
of  the  secondary  signal  at  the  start.  A  digital  computer  solution  will  be 
required  in  most  cases.  The  procedure  is  as  follows: 
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Figure  3.6.  Sketch  of  the  Secondary  Signal  a  (/it)  = 

and  the  Corresponding  Modified  Nonlinear  Element 
for  the  Relay  with  Dead  Band 


Step  1 

Assume  an  approximate  secondary  signal  cr(/3t)  of  the  staircase  type 
shown  in  Figure  3.7.  The  signal  <r(/3t)  has  the  same  properties  which  have 
been  assumed  above  and  described  by  Equation  (3.4) .  The  amplitude  B  is 
divided  into  n  equal  increments  (B  =  nA) .  The  width  of  the  ith  increment  is 


T.  and  is  unknown.  Define  the  T.  so  that 


2A 

A 


r 

s 


.2  Ti  =  */2  •  (3. 1 

Step  2 

Assume  the  composite  signal  Aq  *  as  an  input  to  the  nonlinear 

element  where,  tor  convenience,  the  bias  Ao  is  confined  to  the  A  increments 
as  well.  That  is,  let 


where 


A 

o 


jA  -  B  , 


(3.50) 


J  ~  11  ■  (3.51) 

Both  j  and  n  are  positive  integers.  Now,  write  out  the  algebraic  equation  tor 

the  average  output  of  the  nonlinear  element  over  a  complete  cycle  ft  =  2ir. 

The  equations  take  the  following  form: 
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Av  =  —  |n(Aq  +  A^Tj  +  n(Aq  +  2a)t2  +  ...  +  n(Aq  +  iA^T.  +  ... 

+  N(Ao  ♦  ni)Tn]  +^-[n(Ao  -  a)t,  +  N(Ao  -  2a)t,  ♦  ... 

+  N^Aq  -  kA^Tk  +  ...  +  N^Aq  -  nA^J  ,  (3.52) 


or 


A 

v 


1  n 

ia\t.  +  —Y.  n(A 

/i  !r.“  \  o 


-  iA|T. 


(3.53) 


Substituting  the  value  for  Aq  given  in  Equation  (3. 50)  and  recognizing  that  A^ 

a 

is  really  N  at  this  particular  value  Aq  gives 


N.  =  —  Y  [N(jA  +  iA)  +  N(jA  -  iA)  ] T.  .  (3.54) 

J  *  i-1  1 

a 

N\  is  defined  by  the  equation 

N.  i  N(A  =  jA^  .  (3.55) 

3  \  o  / 

Step  3 

Place  the  desired  values  of  N*  on  the  left-hand  side  of  Equation  (3.54) 
and  solve  for  the  unknown  T..  There  will,  in  general,  be  (n  -  1)  linearly 
independent  equations  and  the  relation  given  in  Equation  (3.49)  is  required  to 
complete  the  solution.  If  no  negative  or  imaginary  values  are  found  for  the  T. , 
this  means  that  secondary  signal "a(/3t)  has  been  found  which  approximates  the 
desired  signal  a*  (/3t) .  Since  the  amplitude  B  of  the  secondary  signal  is 
chosen  a  priori,  it  is  advisable  to  have  a  variable  kj  multiplying  the  desired 

A 

N.  For  example, 

N*  =  kj  IaJ™  Isgn  ^A  \  .  (3.56) 
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f 


k 


i 


4 


To  specify  k4  initially  would,  in  many  cases,  eliminate  any  possibility  of 
obtaining  a  physically  realizable  answer.  It  is  remarked  that  after  the  wave- 

A 

form  of  cr(/3t)  and  the  magnitude  of  N  have  been  determined,  a  second  choice  of 
B  might  lead  to  a  better  approximate  solution. 

Step  4 

Whena(^t)  has  been  found,  draw  a  smooth  curve  through  the  midpoint 
of  each  increment.  This  should  give  the  best  approximate  solution  a*  (  fit)  . 

The  value  of  n  chosen  will  depend  upon  the  desired  accuracy.  It  is  suggested 
that  the  first  try  be  made  with  a  small  n. 

A  simple  example  of  the  linearization  of  the  absquare  will  be  given  to 
show  how  the  method  may  be  used. 

Example  3  (Absquare) 

The  absquare  is  described  by  the  equation 

N(e)  =  e2  sgn  (e)  .  (3.57) 

The  desired  modified  nonlinear  element  is  given  by 

N*|Ao;  a(0tjj  =  kjAQ  ;  Aq  <  B  .  (3.58) 

Since  N(e)  is  described  b>  a  one-term  polynomial,  it  is  not  necessary  to 
choose  a  specific  value  for  B.  Choose  n  =  5  for  this  example  and  the  value  of 
A  becomes  B/5.  Equation  (3.65)  is  written  out  as  follows: 

Ak,  =  (4A2Tj  +  8A2T2  +  12A2T3  +  16A2T4  +  20A2T5)  (3.59) 

2Akj  =---  (9A2Tj  +  16A2T2  +  24A2T3  +  32A2T4  +  40A2T5)  (3.60) 

3Akt  =  (20A2T4  +  40A2T2  +  50A2T3  +  64A2T4  +  60A2T5)  (3.61) 
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1 


(3.62) 


4Ak4  =  —  (34A2T4  +  40A2T2  +  50A2T3  +  64A2T4  +  80A2T5) 


5Akj  =-i(50A2T!  +  58A2T2  +  68A2T3  +  82A2T4  +  100A2T5)  .  (3.63) 


Successive  elimination  of  the  unknown  terms  shows  that 


Tj  =  T2  =  T3  =  T4  -  0 


This  leaves  finally 


(3.64) 


,  20  2 
Ak<  = - A2T5 

1  7T  s 


(3.65) 


From  Equation  (3.49)  it  is  seen  that 

T5  =  tt/2 

Thus, 


(3.66) 


k4  =  2B*  ;  Aq  <  B*  .  (3.67) 

Equation  (3.66)  relates  that  cr*  (j3t)  is  a  square  wave  and  Equation  (3.67) 
states  that  the  resulting  gain  is  2B*  in  the  linear  region.  This  answer  has 
already  been  shown  in  Example  2  to  be  correct.  When  the  nonlinear  element 
is  piecewise  linear ,  the  amplitude  B  does  not  have  to  be  chosen  first  and  the 
parameter  A  may  be  factored  out  of  Equation  (3.54)  in  a  manner  similar  to 
the  above  example. 

til 

The  examples  given  using  the  method  of  specifying  the  m  ‘  derivative 
of  nJa^;  cr(/3t)J  shows  that  only  very  limited  success  can  be  expected  from  any 
method  where  the  nonlinear  element  is  more  general.  This  applies  to  the 
approximate  method  described  above  where  the  nonlinear  element  may  exist 
only  in  graphical  form.  Although  emphasis  has  been  placed  on  the  region 
Aq  ^  B  around  the  origin,  's  restriction  is  not  necessary.  The  region  of 


\ 


» 
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N  for  Aq  >  B  may  sometimes  be  altered  in  a  specified  manner  by  the 
presence  of  a  (ftt) . 

If  the  desired  secondary  signal  cannot  be  found  by  any  of  the  methods  des¬ 
cribed  above,  there  is  an  alternative  method  of  modifying  the  nonlinear  element 
which  may  prove  effective  [suggested  to  the  author  by  C.  D.  Johnson].  Inorderto 
apply  this  alternative  method,  the  desired  K  *  or  N*  must  lie  in  the  region  R  which 
is  attainable  with  some  secondary  signal,  it  is  usually  not  difficult  to  establish  R. 
The  first  step  is  to  choose  a  convenient  secondary  signaHH  (h) ,  perhaps  a  sine  wave. 
Step  2  is  to  determine  how  the  amplitude  B  of  ai  Pt)  mustbe  varied  in  accordance 
with  the  input  signal  amplitude  Aq  in  order  to  give  the  desired  modified 
nonlinear  element.  This  may  be  done  by  solving  the  equation 

Nj/-o;  cr(/?t)J  =  N*(A)  .  (3.68) 

The  unknown  in  Equation  13.68)  is  B  and  will  only  be  function  of  A^  once  a(pt) 
is  chosen.  This  function  B(A)  will  then  become  a  "variable  gain  box"  which 
will  premultiply  cr(pt)  before  it  is  summed  with  the  input  A^.  From  a  hard¬ 
ware  implementation  viewpoint,  this  alternative  method  is  not  appreciably 
more  complex  than  shaping  the  secondary  signal  o*(  pt)  in  the  inverse  DJDF 
method. 

Oldenburger  and  Ikebe  [40]  have  offered  still  another  linearizing 
method  where  a  relay  function  and  a  triangle  wave  secondary  signal  are  placed 
in  series  ahead  of  the  existing  nonlinear  element.  In  effect  the  insertion  of 
the  relay  function  in  front  of  another  nonlinear  element  in  the  signal  flow  path 
results  in  a  relay  function  for  the  combination.  It  was  seen  in  this  chapter 
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that  the  relay  is  perhaps  the  most  versatile  nonlinear  element  when  one  is 
concerned  with  shaping  the  modified  nonlinear  element  with  an  extra  signal. 
Therefore,  the  approach  presented  in  this  chapter  can  be  used  to  extend  the 
method  of  Oldenburger  and  Lee  be  to  the  more  general  problem  of  specifying 
nJaq;  (r(£?t)j rather  than  strict  linearization. 


C  HA  PI  ER  IV 


THE  MODIFIED  NONLINEAR  ELEMENTS  AND  DUAL  INPUT  DF'S 
FOR  TEN  NONLINEAR  ELEMENTS 

In  this  chapter  the  modified  nonlinear  element  and  the  DIDF's  (Chapter  II) 
of  ten  common  nonlinear  elements  are  given  for  the  case  when  the  secondary 
signal  is  either  a  sine  wave,  a  triangle  wave ,  or  a  square  wave.  These  secondary 
signal  waveforms  were  chosen  to  show  the  qualitative  trend  in  N[A;o-'8t)  ]  as 
the  waveform  is  changed.  The  specific  nonlinear  elements  are  as  follows: 

a.  Absquare 

b.  Relay 

c.  Relay  with  dead  band 

d.  Preload 

e.  Dead  band 

f.  Limiter 

g.  Limiter  with  dead  band 

h.  Relay  with  dead  band  and  hysteresis 

i.  Relay  with  hysteresis 

j.  Limiter  with  hysteresis. 

Closed  form  solutions  for  the  modified  nonlinear  elements  are  given  for  all 
cases.  Closed  form  solutions  are  also  given  for  the  DIDF's  with  the  exception 
of  five  nonlinear  elements  with  sine  wave  secondary  signals.  Approximate 


mm 


solutions  are  given  for  some  of  these  cases  in  the  form  of  finite  series 
obtained  by  trapezoidal  integration  tabulated  along  with  the  complete  set  of 
curves  of  both  the  modified  nonlinear  elements  and  the  DIDF's.  The  average 
DIDF  method  as  explained  in  Chapter  II  was  used  to  calculate  the  TSIDF  with 

the  exception  of  the  three  memory  type  nonlinear  elements.  Otherwise,  the 
modified  nonlinear  element  method,  also  explained  in  Chapter  II,  was  used  to 

calculate  the  DIDF.  When  the  modified  nonlinear  element  is  derived,  no 
special  effort  is  made  to  insure  that  the  equation  holds  for  negative  inputs  A  . 
It  is  simply  necessary  to  recall  that  the  modified  nonlinear  element  is  an  odd 
function. 


It  is  emphasized  again  that  the  DIDF's  of  the  double  valued  nonlinear 
elements  found  by  the  modified  nonlinear  element  method  are  valid  only  when 
the  frequency  of  the  secondary  signal  is  very  high  compared  to  that  of  the 
primary  input  component.  An  example  is  given  in  the  appendix  to  show  the  pro¬ 
cedure  used  in  deriving  these  modified  nonlinear  elements  and  DIDF's. 
Numerical  answers  for  the  DIDF's  were  obtained  by  digital  computer  evalua¬ 
tion  of  the  given  equations. 

The  following  definitions  describe  angles  and  symbols  used  in  the 
modified  nonlinear  element  and  DIDF  equations: 

IXl  —  absolute  value  of  X  (4.1 


sgn  (X) 


1;  X  >  1 


-1;  X  <  1 


sat  (X) 


X;  IXl  £  1 


sgn  (X) ;  IXl  >  1 


* 
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4.1.  Absquare 


N(e)  =  ez  sgn  (e) 


(4.  5] 


Figure.  4.1.  Absquare 

Sine  Wave  Secondary  Signal  (Figures  4.2  and  4.3) 


5(vb)  =  ' 


(“o  +  b!)s1“J(4)  +  3Ao BiJ  -  (if)  .  •  Ao 
l(V+H  >  Ao£B 


K(A,B)  =  [(7AZ  +  B<:)E(k)  -  (A  -  B)2F(k)]  ,  (4.7 


=  SJA+J}  f^7A2  +  B2) 
2\TaB 


<  B 


(4.6; 


k  = 


(A  +  B) 

Triangle  Wave  Secondary  Signal  (Figures  4.4  ar.d  4.  5) 


A  3 

A  B  +  ~  ;  A  £  B 
o  .7B  o 


(4.8; 


N  A  ,  B  =  \ 


2  ■ 


\o  +  T  5  Ao  >  B 


(4.9; 


2F(k)  denotes  the  complete  elliptic  integral  of  the  first  kind  and  E(k) 
denotes  the  complete  elliptic  integral  of  the  second  kind.  The  modulus  is 
defined  as  k. 


ve  Sec 


Figure  4.4.  Modified  Nonlinear  Element  for  Absquare,  Triangle  Wave  Secondary  Signal 


I 


I! 


K(A,B) 


(t  [(■»  *-£}■"-  ©  •  fl(£  .?*)]. « 


B+W;  A  s  B 


>B 


(4.10) 


Square  Wave  Secondary  Signal  (Figures  4. 7  and  4.  8) 
S(A„'B)  - 

f[.  i(».2)jr© 


2BA  ;  A  <  B 
o  o 


A  +  BJ  ;  A  >  B 
o  o 


(4.11) 


Kw.B)  - 


A  >  B 


2B  ;  A  ^  B 


(4.12) 


4.2  Relay 

N(e)  =  M  sgn  (e) 

Sine  Wave  Secondary  Signal  (Figures  4.  9  and  4. 10) 

I?(Ao’B)=T!si"'‘sat(4) 


(4.13) 


4M  1 


K(A, B)  =  ^2  ^2  t(A  -  B)K(k)  +  (A  +  B)E(k)]  ; 


(4.14) 

(4.15) 


Figure  4. 6.  Relay 


N 

M 

• 

•M 
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Figure  4.7.  Modified  Nonlinear  Element  for  Absquare,  Square  Wave  Secondary  Signal 


‘jrtfp'i’iir  • 


where 


! 


k2  = 


4AB 


(A  +  B) 

Triangle  Wave  Secondary  Signal  ( Figures  4. 12  and  4. 13) 


=  Msat(-r) 


-  M  4M 

K(A,B)  =  —  (2y,  -  sin  2y,)  +  —  cos  7, 


Square  Wave  Secondary  Signal  (Figures  4. 14  and  4. 15) 


(  0  ;  A  <  B 
o 


M  ;  A  a  B 
o 


*(v0 = 

' It I1  -  sm(a)‘ 


4.3  Preload 


N(e)  =  Ce  +  M  sgn  (e) 


(4.16) 


(4.17) 


(4.18) 


(4.19) 


(4.20) 


(4.21) 


Figure  4. 11.  Preload 
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N(A  ,  B) 


Figure  4. 12.  Modified  Nonlinear  Element  for  Perfect  Relay,  Triangle  Wave  Secondary  Signal 


A 


Figure  4.13.  DIDF  for  Relay,  Triangle  Wave  Secondary  Signal 
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.ivV&i 


C,M) 


Figure  4. 16.  Modified  Nonlinear  Element  for  Preload,  Sine  Wave  Secondary  Signal 


18.  Modified  Nonlinear  Element  for  Preload,  Triangle  Wave  Secondary  Signal 


AC 

M 


5  4.19.  DIDF  for  Preload,  Triangle  Wave  Secondary  Signal 
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C,  M) 


for  Preload,  Square  Wave  Secondary  Signal 


1 


I 


f 


l 

c 

I 

I 


K(B,B,a)  =  -  -  M 
tt2B 


where 


(A  +  B)2  -  a2 

(A  +  B)2 

r 

(4.33) 

LE(4,k)  "  4B2  F(l’k)J 

;  A  =  B  (4.34) 

k2  =  i  _  jL 

4B2 

(4.35) 

I 

1 

:  Sin-'  r  f  , 

*J4B‘  -  a2 

i 

1 

(4.36) 

Case  c 


4M 

K(A,B,a)  =  -J~2  - 

*A 


(A  -  B) 2  > 


2  .  „2 


+  B)2  -  a2 
2 


(A2  -  B2) F  +  [(A  +  B)2  -  a2]E 


(2A  ~  a)  (^T-  ljnfa2,k) 


(4.37) 


where 


2  4AB 
a  =  - 


®i2  = 


(A  +  B)2 
2 


A  fa2  -  A2  -  B2  \ 

B  \  2A2  -  a2  / 

Triangle  Wave  Secondary  Signal  ( Figures  4 .  25  and  4 .26 ) 


(4.38) 

(4.39) 


T  /a  +  A_ 


a  -  A 


-  sat 


”(Ao'B'a)  =t[ 

A  M  f 

K(A,  B,a)  =  — |y3  -  y2  -  sin  y,  cos  y3  +  siny2cosys 

.  ~  B)  4g  "I 

A  (cos  y3  -  cos  y2)  +  — cos  y3 


(4.40) 


(4.41) 
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Figure  4.25.  Modified  Nonlinear  Element  for  Relay  with  Dead  Band,  Triangle  Wave  Secondary  Signal 


I 


gure  4.  28.  Modified  Nonlinear  Element  for  Relay  with  Dead  Band,  Square  Wave  Secondary  Signal 


*V| 


Modified  Nonlinear  Element  for  Dead  Band,  Sine  Wave  Secondary  Signal 


f 


l  i 


»  K 


* 


where 


«i‘  (1*loos*i)  -^n-'  sat  (-57) 
+  Csati  -Is 

V  D.  D. 

'■»i 


1  -  sat  ( — 


1 '  *1 


(i) 


(4.47) 


and 


D.  ^  ^A2  +  B2  -  2AB  cos  <+>. 

/a2  -  A2  -  R2) 
V  2AB  > 


$  ~  cos-1  sat 

s 


i  ,  i  -1  4  /a2  -  A2  -  B2\ 

d>.  =  —  rf>  = —  cos  sat  l - — — - / 

vi  n  ^s  n  \  2AB  ' 


n  =  integer  >  2  . 

The  accuracy  is  increased  as  n  is  allowed  to  become  larger. 

Triangle  Wave  Secondary  Signal  (Figures  4.32  and  4.33) 

N(A  ,B,a)  =—  —  (eo 2  -  0o J)  +  A  (*  -  0n-  -  60J\ 
V  o  /  7r  7r  \  3T  2T  /  o\  2T  3T ) 

T  a(02T  “  63t) 


(4.48) 

(4.49) 

(4.50) 

(4.51) 


n(a  ,B,a^  =  CA  fl  +  ^  ;  B  >  A 

\  o  /  0  \  B  /  o 


(4.52) 

(4.53) 


*  C  R 

K  (A ,  B,  a)  =  —  T 
7T  A 


-  *)*  *  (4)! 


(cosy2  -  cos  y3) 


+  sin  y3  cos  y3 


/A  aA\ 
\B  +  B2" ) 


+  sin  y2  cos  y2 


(A  aA  \ 
\B  B2 ) 


_4a_ 

B 


cos 


(A  aA\  (A  aA\ 
Y3  “  Y3  \  B  +  B2  )  ~  Yz  \  B  ”  B2  / 


l/Ay(cosJy3  -  cos3  \2)  +— 7 r 
T \B/  B 


K(A,B,a)  =  C 


B  >  A  +  a 


(4.54) 

(4.55) 
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Modified  Nonlinear  Element  for  Dead  Band,  Triangle  Wave  Secondary 


Figure  4.35.  Modified  Nonlinear  Element  for  Dead  Band,  Square  Wave  Secondary  Signal 


f 


gure  4.37.  Modified  Nonlinear  Element  for  Limiter,  Sine  Wave  Secondary  Signal 
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Square  Wave  Secondary  Signal  (Figures  4.42  and  4.43) 


N^A^.B.a^  =-^-|b^2  +  sgn  ^a  -  B  -  -  sgn  ^a  +  B  -  A^  j 

+  ^Aq  -  a^  |"sgn  ^a  -  B  -  +  sgn  ^a  +  B  -  A^j  J 

+  ^Aq  +  a  -  B^  j^l  -  sgn  ^  B  -  a  -  Aq 


(4.69) 


K(A,  B,  a)  =  — 


y3  +  y?  -  sin  7s  cos  Y;  -  sm  y2  cos  y2 


2B  2a  1 

+  —  (cos  y3  -  cos  y2)  +  —  (cosy3  f  cos  y-.)J  (4.70) 


4.7  Limiter  with  Dead  Band 

N(e)  =  cjjbsat^-^  -  asat^^J  (4.71) 


Figure  4.41.  Saturation  with  Dead  Band 


Sine  Wave  Secondary  Signal  (Figures  4.44  and  4.45) 

N(Ao>B,a,b)  =—  [((?*  +  ~  n2  "  03) Aq 

+  B(cos  05  -  cos  Oq  -  cos  03  +  cos  02) 
+  b(#5  -  @,j)  -  a (03  -  (?2)  j 


(4.72) 


I 


"  Limiter,  Square  Wave  Secondary  Signal 


Modified  Nonlinear  Element  for  Limiter  with  Dead  Band,  Sine  Wave  Secondary  Signal 


K(A,B,a,b)  =  ~f  J^1  “  +  (4)  ]  (cos^2  -  cos  73) 

'  [i1  '!)  +  (I)  ]  (C0S74  “  C0S^) 

/A  aA\  .  /A  bA\ 

+  sin  73  cos  y3  "  sin  75  cos  (5  (-  + 

/A  aA\  .  /A  bA\ 

+  sin  72  cos  72  ^  -  ^7 )  ~  sln  74  cos  74  -  -gj ) 


4b  / 

A  +  aA 

\ 

(A 

—  cos  75  -  73^ 

B  B7 

)  + 

Hi 

/A  aA' 

1  / 

A_ 

bA\ 

72  \  B  ”  "b7, 

)  +74( 

B  “ 

i7) 

1  /A\2 

+  ll\B/  (cos3-y3  "  cos3  72  -  cos3  75  +  cos3  74)  (4.84) 

Square  Wave  Secondary  Signal  \  Figures  4.48  and  4.49) 

n(Aq,  B,a,b)  ■=  sgn  (b  -  B  -  A^-  sgn  (b  +  B  -  A 

-  sgn  (a  -  B  -  A^j  +  sgn  (a  +  B  -  A^ 

+  (  i-2)  [sgn  (a  -  B  -  Ao)  +  sgn  (a  +  B  -  Ao) 

-  (^r2)[sgn  (b  -  B  -  Ao)  +  sgn(b  +  B  -  Ao)J 

+  ( - — °)  [l  -  sgn  (B  -  >  -  Ao)J 

-  ( - b— 2)  [l  -  (B  -  b  -  A0)  1  1  (4.85) 


K(A,B,a,b)  = — y5  +  74  “  73  -  72  +  siny2  cos  72  +  sin 73  cos  73 


-  sin  74  cos  74  -  sin  75  cos  75  +  ——  (cos  72  -  cos  73 


-  cos  y4  +  cos  75)  -  —  (cos  72  +  cos  73) 
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Figure  4.49.  DIDF  for  Limite: 

Square  Wave  Seconda: 


f 


2b  ,  .1 

+  —  (COS  74  +  COS  7S) 


4.  8  Relay  with  Dead  Band  and  Hysteresis 


(4.86) 


H(e) 

LA 

-b  -c 

m 

0  1 

i  e 

— 

-M 

Figure  4.50.  Relay  with  Dead  Band  and  Hysteresis 

Sine  Wave  Secondary  Signal  (Figures  4.51  and  4.52) 

*  M  v-  _  o 

N(A,B,a,b)  =  —  (0S  +  03  -  04  -  42)  ;  B  >  ^ —  ,  (4.  87) 

subject  to  the  restrictions  that 

9  2  is  set  to  — 7r/2  when  94  -  -tt/2 
04  is  set  to  +7r/2  when  02  =  tt/2 
03  is  ."-et  to  +7t/2  when  05  =  ir/2 

05  is  set  to  -tt/2  when  fl3  =  -it/ 2  (4.83) 


-  2m  "%4  r  ,  / 

K(A,  B,a,b)  =  -j-  J  sin-1 


.  t  /a  +  A  sin  cut \  .  _i 

sin  I - — - 1  +  sm 

(h  +  A  sin  cut\ 

V  B  / 

\  B  / 

sin  cut  dcut 


y3r 

2M  r  .  ..J 

~7aJ  sin  ( 

n 


a  -  A  sin  cut\  .  _<  (  b  -  A  sin 
- - - )  +  sin  l - - - 


in  cut\  . 

- jjsi 


sin  cut  dcut 


*xvl 

—  (COS  y4  +  COS  y3) 


(4. 89) 
ill 


N(*> 


These  integrals  may  be  solved  in  closed  form.  There  are  two  types  to  con¬ 
sider.  Replacing  a  and  b  by  c  to  give  more  generality,  the  two  integral  fypes 
become 


2M~f*  .  I 

K>  ’7Z  J  sm  V 


c  +  A  sin  cot 
B 


sin  cot  dcot 


(4.90) 


K, 


2M  Jr  .  / 

=  7iJsin  V 


c  -  A  sin  cot\  . 

— - - )  sin  wt  do)t 

B 


(4.91) 


where 


•rr  f_l  .  1  /C\  cos  y4  _!  (c  -  A  sin  y4\ 

K'  =  S  [lSln‘  (b)  -  -X^Sm  { - B - -)_ 

K,  =  ™[f  sin-1  ^  sm-1 


+  K 


li 

14. 92) 


K2j 


(4.93) 


Yl 

Kll  =  / 
o 


COS2  cot  dcot 


1  - 


c  +  A  sin  cot\ 2 
B 


(4.94) 


r  2 


K21  =  f 


cos2  cot  dcot 


-  A  sin  cot\ 2 


(4.95) 


B 


A  general  form  of  the  solution  for  both  K41  and  K2j  is  given  as  follows: 


K  = 


B. 

Thg 


^y[E(^.k)  +  F(0,  k) 


x  2(k2 

/  9  k2  \  n 

+  (a2  -  2  +~'n(<£,a2,k) 


a2  sin  cj>  cos  <j>  4  1  -  k2  sin2  <p 

1  -  a2  sin2  <p  J 


(4.96) 

(4.96) 


The  correct  values  of  k,  g,  a2,  <p,  and  h  are  given  in  Table  4-1. 
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\ 


N(A,B,a,b)  = 


Relay  with  dead  band  2  ( a  +  b) 
and  hysteresis  (b  -  a  -  2B) 
(Figure  4. 53) 


*4.97) 


for  B  < 


b  -  a 


Figure  4.53.  Modified  Nonlinear  Element 


lK(A,B,a,b)  !  =-^WlPr 

7T  A 


(Tj  +  Y4/  .  B  <  (4.98) 


=  smy,\  (y.-Ti),  p  .  b  -  a 

—  \cosyi  +  cos  y2  /  2  ’  2 


(4.99) 


Triangle  Wave  Secondary  Signal  (Figures  4.  54  and  4.55) 


Case  a 


Both  N 


(Ao-B>a.b) 


and  K(A.. B,a,b)  are  the  same  as  in  the  sine  wave  case. 


Case  b 


b  -  a 


<  B  <  b 
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(4.105) 


K(A,B,a,b)  =  (same  as  case  b) 
Square  Wave  Secondary  Signal  (Figures  4.57  and  4.58) 
Case  a 

/B  <  b  -  a  \ 

\B  < 

See  Case  a  for  sine  wave  secondary  signal. 

Case  b 

(> 


B  > 


b  - 


a  r 

N(Ao’B,a,b)  =TLSgn(Ao  +  B  +  a)  +  sgn(AQ  +  B  -  b) 


sgn 


-  B  +  b)  +  sgn(AQ  -  B  -  a)J 

(4.106) 

2M  r 

-  —  [cos  yi  +  cos  ys] 

(4.107) 

Figure  4.  59  shows  the  phase  of  DIDF  for  relay  with  hysteresis  and 
dead  band. 


4.  9  Relay  with  Hysteresis 


Figure  4.56.  Rdlay  with  Hysteresis 


Square  Wave  Secondary  Signal 


ith  Dead  Band  and  Hysteresis, 
condary  Signal 


Sine  Wave  Secondary  Signal  (Figures  4.60  and  4.61) 


H(.) 


Figure  4.60.  Modified  Nonlinear  Element  for  Relay  with  Hysteresis,  Sine  Wave  Secondary  Signal 


I  > 


Figure  4,  64.  Modified  Nonlinear  Element  for  Relay  with  Hysteresis,  Square  Wave  Secondary  Signal 


(fe  Sr  g&ft-'i 


Figure  4.69.  Phase  of  DIDF  for  Limiter  with  Hysteresis, 
Siue  Wave  Secondary  Signal 


igure  4.  r,2.  Modified  Nonlinear  Element  When  B 


N5  =  j  A  (7T  +  2 02)  +  7Tb  -  2a 0o  - 


2  -  2B  cos  02 


B  -  s  A  . 

o 


K(A,  B,a, b)  jBJl  -  (bl.Asino,.y4B  ^  A^sin  ^ 

+  A  sin  ujt  |gin-i^a  ~  A^sin  R1-n-i^b  ~  A  sin  cjtj 

+  sin-Y-l-  ABSin  “l)  +  o.„_^b  +  A  sin_cot\ 

+  b  sin-i  (  bjj^ sin  ut_j  _  b  s.n_t  ) 

+  a  sin-1  f  —  A  sin  —  ^  _  o  /a  -  A  sin  wt\ 

V  B  /  a  sin  y - - - j  -sinwtdwt 

y 3  i 

+  7T  I  Ja  Sin  wt  ("tt  +  2  sin"1  A  sin  \  1 


+  A  sin  cjt 


+  sin" 


+  a  sin" 


a  -  A  sin  cjt\  | 

g - )  /  sin  cot  dcot 


+  rb  -  2a  sin"1  ('»-  -  A  sin  cot  \ 
'  B  / 


-  2B  Jl  -  (— V"  “-)  sin.td.t+iMcos 


ACos^3  •  (4.124) 


Triangle  Wave  Secondary  Signal  (Figures  4.73  and  4.  74) 
Case  a  (Figure  4.71) 


where 


(B  <  b) 


T^1  +l)A0  -  T-  -  8  +  2b  -•£]  ; 


a  "  B  £  A  Sa+B 
o 


K(A, B, a, b)  =  /g  2  +  h  2 
V  A  A 


K  =  tan-1 — 


=_c_r 

gA  ttA  [A(y4  +  T2  -  sin  y2  cos  y2  -  sin 


74  cos  y4 


+  2b(cos  y3  +  cosy4)  +  a(cos  y2  +  cosy.) 


(4.125) 


(4.126) 


(4.126) 


(4.127) 


A  /  a\ 

2  \  b  /  'yz  +  y3  +  sin  y2  cos  y2  -  sin  y3  cos  y3) 

A2  / 

2B  (cos73  -  COS  y2  +_cos3y3  -i-COS3y2) 

7  (B  +  i* )  (cos  73  -  cos  y2)  ;  B<-^-i  (4>128) 


Modified  Nonlinear  Element  for  Limiter  with  Hysteresis,  Triangle  Wave  Secondary  Signal 


and 


7rA  ^  (sin2  y  \  -  sin2  y2)  ~  2bsiny4  +  a(siny2  +  sin  y3) 

+  T  i1  +  £){sin2V2  -  sin2Y3)  (sin3y2  -  sin3  y3) 

"T  (B  +  B~)  (siny2  -  slnV3)  ;  B  <  -  ~  a  ,  (4.129) 

C  f 

»A  “  ttA  lA  (y2  +  y4  -  sin  y2  cos  y2  -  sin  y4  cos  y4  +  2b  cos  y4 

+  a(cosy2  +  2  cos  y3  -  cos  y4)  +  A  (l  +  ^)(2y3  -  y4 

-  y2  +  sin  y2  cos  y2  -  2  sin  y3  cos  y3  +  sin  y4  cos  y4) 

A2  /  i  i 

“  ~2B  icos  72  “  2c°sy3  +  cos  y4  cos3y2  -_cos3 

+  — cos3y3)  +  ^B  + -g*)  (2  cos  y3  -  cos  y4  -  cos  y2) 


y  i 


;  (4.130) 


b  -  a 


£  B  s  b  , 


^A  ttA  74  ~  sin2y2)  -  2bsiny4  +  a(siny4  +  siny2) 

+  _2~(1  +  B)(sin2y2  “  sin2  74)  "  ^  (sin3  y2  -  sin3y4) 

1  /  a2\ 

'TV6  +  B-J<Siny2  "  Sinv4) 


b  -  a 


B  <  b  .  (4.131) 


Case  b  (Figure  4.62) 


(B  >  b) 
A 


N, 


A  <  B  -  b 
o 


(4.132) 
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^=T\?V+b)A0 


-  B  +  2b  -  —  J  ; 


(4.133) 


B-b^A  ^  B  +  a 
o 

N  (a  ,B,a,b,c^  =  M  ;  A  >  B  +  a 

'0  /  o 

K(A,B,a,b,c)  =—  ( — g—)  (^'4  +  sin y4  cos  y4) 

+  (1  +  +  74  "  sin  y4  cos  y4 

A  / 

-  sin  73  cos -y3)  -  —  ^cos  y4  -  cos  y3 

1  3  1  3  \ 

+  Ycos  73  --3  COS3  74) 

1  /  a2\ 

+  —  (2b  -  B  -  -  )  (cos  74  -  cos  73) 

2  1 

+  "A  (b  +  C0S  B  >  b 

Square  Wave  Secondary  Signal  (Figures  4.75  and  4.76) 


(4.134) 


(4.135) 


Case  a  ( Figure  4.77) 


(B  <  b) 


N3  =-f-(A0  +  b  -  b)  ;  B-b<AQ<B  +  a  (4.136) 

c  Ta 

"  ?rA  7  (2yi  +  72  +  “  2  sin  74  cos  74  -  sin  y2  cos  y2  -  sin  y3  cos  73) 

-  2b  cos  73  +  a  (cos  y2  +  cos  73)  +  B(cosy3  -  cosy2)l  ; 


Figure  4.  75.  Modified  Nonlinear  Element  for  Limiter  with  Hysteresis,  Square  Wave  Secondary  Signal 


+  B  b  .  B 


a  -  B '  |  b  -  o 

B<  -r— 


I  II2 


a  /b  -  B  a  +  B  A 

f  i  ; 

1  1  b  -  o 

I  I  —  'B'k 

I  I  I 


Figure  4.77.  Modified  Nonlinear  Element  When  B  <  b 
for  the  Square  Wave  Secondary  Signal 


C  ["a 

hA  =lrk\2(2  sin274  "  sin2'>'2  '  sift2  y3)  +  b(siny4  -  siny2) 
+  a  (sin  y4  +  siny3)  +  B(siny3  -  siny2)J  ; 


b  -  a 


(4.138) 


C  [A 

-  ^  r^('y4  +  72  +  2y3  -  sin  y4  cos  y4  -  sin  y2  cos  y2  -  2  sin  y3  cos  y3) 


+  — (3  cos  y4  -  cos  y2  -  2  cos  y3)  +  a(cos  y2  -  cos  y.) 

2  -J 


b  -  a 


<  B  <  b  , 


(4.139) 


C_[A 

ttA  (_  2 


(sin2  y4  -  sin2  y2)  +— (3  siny2  -  2  siny3  -  3  siny4) 


<X 

— (3siny2  -  2  sin  y3  +  siny4)  +  B(siny4  -  siny2) 


b  -  a 


<  B  <  b  . 


<4. 14(B 


Case  b  (Figure  4.78) 


(B  >  b) 


Figure  4.78.  Modified  Nonlinear  Element  When  B  >  b 
for  the  Square  Wave  Secondary  Signal 

N(A,  B,a,b)  (B  +  a)sat  (— °  J  -  (B  -  b)  sat  (-  0  ) 

'  B  b }  __ 

(4.1 

9C  A 

K(A ,  B,  a,  b)  -  —  ~(y3  +  y4  -  sinv3  cos  y3  -  siny4  cos  y4) 
b  cos  y4  +  a  cos  y3  +  B(cosy3  -  cosy4)l  ; 


B  >  b 


(4.1 4 
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CHAPTER  V 


TWO-SINUSOID  INPUT  DF  FOR  RELAY  WITH  HYSTERESIS 


In  this  chapter  the  TSIDF  is  calculated  for  the  relay  with  hysteresis, 
taking  into  account  the  frequency  ratio  of  the  input  components.  This  is  done 
by  using  the  "Average  DIDF"  method  explained  in  Chapter  n.  The  nonlinear 
element  is  defined  in  Figure  5. 1  and  the  essential  mathematical  steps  are 
indicated  in  the  following  equations: 


1 

1 

N(e), 

1 

M 

i 

®Q 

1 

1 

0  e 

•M 

Figure  5.1.  Relay  with  Hysteresis 


K  -  K  (A,  B,  n) 
r  r 


2  2n  2n 

717  /  d<£  /  N[A  sin  pwt  +  B  sin  (qwt  +  d>)  ]sin  pwt  dwt 
2ff  A  J  J  ,,  .. 

o  (5.1) 
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r 


2tt  2tt 


K.  =  K.(A,B,n)  =  T”27  f  dip  f  N[A  sin  pcut  +  B  sin  (qu;t  +  <p)  ]cos  pcut  doit 
ll  2tt  A J  J  _ 

o  o  (5.2) 


K 


!  =  Ik  2  +  k.2 

v  r  l 


K  =  tan"1  — 


K. 

_ l 

K 


(5.3) 


(5.4) 


where  K  ,  K.,  Kl,  and  /K  are  the  real  part,  imaginary  part,  magnitude,  and 
phase,  respectively,  of  the  DIDF.  The  primary  component  in  the  input  signal 
is 

ej  =  A  sin  poit  (5. 5) 

aiid  the  secondary  component  of  the  input  signal  is 

e2  =  B  sin  (qwt  +  <p)  (5.6) 

It  is  not  necessary  to  identify  input  frequencies  p  and  q  since  only  the  fre¬ 
quency  ratio  n  =  q/p  is  important.  However,  they  have  been  included  in  the 
formulation  for  computational  purposes.  This  will  be  explained  further  in  the 
discussion  which  follows. 

It  was  stated  in  Chapter  II  that  the  average  DIDF  as  computed  for 
multivalued  nonlinear  elements  has  meaning  only  when  n  is  an  irrational 
number.  However,  only  integers  are  chosen  for  p  and  q  in  this  analysis  and 
it  must  be  kept  in  mind  that  the  resulting  DIDF  holds  only  for  irrational  fre¬ 
quency  ratios  in  close  proximity  to  the  chosen  n.  Although  a  closed  form 
solution  can  be  found  for  the  DIDF  when  n  equals  one,  it  is  doubtful  that  the 
integrals  of  Equations  (5. 1)  and  (5.2)  can  be  solved  analytically  when  n  is 
other  than  one.  It  would  be  an  important  contribution  if  only  the  first 
indicated  integration  could  be  performed  analytically  as  in  the  modified 


*. 
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nonlinear  element  method.  However,  it  appears  that  even  this  cannot  be  done 


and  for  this  reason  a  double  integration  by  some  numerical  technique  is  usually 
required.  Therefore,  when  a  DIDF  is  calculated  with  a  digital  computer,  one 
can  expect  to  use  a  relatively  large  amount  of  computer  time. 

In  this  chapter,  the  results  of  a  digital  computer  solution  are  presented 
in  the  form  of  graphs.  Nine  values  of  n  and  four  values  of  B/a  are  chosen. 

A  table  of  seventeen  <p' s  was  set  up  for  "averaging"  the  DIDF  and  performing 
the  second  integration.  The  Awt  increments  used  for  carrying  out  the  first 
integration  were  approximately  one  degree.  It  was  found  that  the  size  of  the 
A<p  increments  did  not  affect  the  accuracy  of  the  final  answer  to  the  same 
extent  as  the  size  of  the  Aojt  increments.  The  increment  for  the  second 
integration  was  chosen  as  A<p  =  27r/16  because  of  limited  available  computer 
time. 

Several  restrictions  must  be  kept  in  mind  when  working  with  the  multi¬ 
valued  nonlinear  elements.  First,  the  nature  of  the  nonlinear  element  requires 
that  each  succeeding  wt.  increment  in  the  integration  routine  be  larger  than  the 
previous  one;  that  is,  wt  must  always  increase  (decrease).  Some  integra¬ 
tion  routines  are  thus  not  suited  for  integrating  multivalued  functions.  Second, 
a  different  answer  is  obtained  for  any  one  particular  integration  depending 
upon  the  assumed  initial  output  state  of  the  nonlinear  element  when  the  integra¬ 
tion  limits  are  0  to  2ir.  This  is  apparent  even  in  sinusoidal  DF  analysis  where 
the  correct  procedure  is  obvious.  The  correct  procedure  in  the  DIDF  case 
required  finding  the  angle  wtj  (0  ^  wt]  ^  27r)  where  the  first  "switch,"  if 
any,  occurred  at  the  nonlinear  element  output  (Figure  5.2).  The  integration 
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N(e) 


interval  in  this  ease  was  cot  =  cot-, ,  to  2n  +  cot,.  There  are  certain  values  of 
A,  IT  and  $  for  which  cot,  cannot  be  found  and  the  integral  is  zero.  The  pro¬ 
gram  was  also  run  with  the  integration  limits  cot  =  0  to  cot  =  2i r,  where  the 
hysteretic  relay  output  was  set  to  -M  at  the  beginning  of  each  integration 
cycle.  Some  discrepancy  was  noted  for  the  parameter  ranges  1  <  A  +  B  <  2 
where  n  was  close  to  one.  Otherwise  both  procedures  gave  essentially  the 
same  results. 

Fourier  Series  analysis  requires  that  the  integration  limits  be  over  an 
entire  period.  Thus  when  n  ^  1 ,  p  was  set  to  one  and  q  was  made  to  take  on 
integer  values.  When  n  <  1,  q  was  set  to  one  and  p  was  made  to  take  on 
integer  values.  The  resulting  DIDF's  are  shown  in  Figures  5.3  to  5. 11.  It  is 
seen  that  the  magnitude  of  the  DIDF  shows  relatively  little  dependence  on  the 
frequency  ratio  n  when  n  >  3.  However,  the  phase  angle  of  the  DIDF  is 
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_A_ 

o 

Figure  5.3.  Magnitude  of  DIDF  of  Relay  with  Hysteresis  When  B/a  =  0.5 
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Phase  of  DIDF  of  Relay  with  Hysteresis  When  B/a 


Figure  5.8.  Magnitude  of  DIDF  of  Relay  with  Hysteresis  When 


fix 


91*0 


Figure  5.11.  Phase  of  DIDF  of  Relay  with  Hysteresis  When  B/ a 


strongly  dependent  upon  n,  especially  when  B/a  <  4.  It  is  noted  that  when  n 
becomes  sufficiently  large,  the  DIDF  of  this  method  approaches  that  obtained 
with  the  modified  nonlinear  element  method.  In  that  case  the  phase  shift 
across  the  nonlinear  element  goes  to  zero  for  the  relay  with  hysteresis  when 
B/a  1.  It  has  been  implied  by  previous  investigations  [10,40,41]  that  the 
modified  nonlinear  element  is  valid  for  multivalued  nonlinear  elements  when 
n  2:  10.  The  DIDF's  given  in  this  chapter  show  that  a  more  severe  restriction 
is  required  in  some  cases.  For  example,  when  B/a  =  1,  A/a  =  1.5,  and 
n  =  10,  the  average  phase  shift  of  the  DIDF  is  -11.4  degrees,  whereas  the 
modified  nonlinear  element  method  predicts  that  the  average  phase  shift  will 
be  zero. 

The  DIDF  presented  in  this  chapter  in  graphical  form  does  not  agree 
with  the  one  given  by  Mahalanabis  and  Oldenburger  [31]  because  the  frequency 
ratio  n  does  not  appear  in  their  DIDF.  Therefore,  the  DIDF  of  Mahalanabis 
and  Oldenburger  is  adequately  presented  by  only  one  graph  sheet  showing  the 
magnitude  and  one  graph  sheet  showing  the  phase.  The  authors  imply  that  the 
frequency  of  the  secondary  signal  may  be  either  higher  or  lower  than  the  pri¬ 
mary  component  frequency.  It  appears,  however,  that  their  approximate  DIDF 
expression  is  only  correct  for  a  band  limited  noise  secondary  signal  with  a  low 
frequency  spectrum.  Their  phase  curves  appear  to  compare  with  those  given 
in  this  chapter  when  n  is  somewhere  between  1  and  2. 
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CHAPTER  VI 


CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 


The  author's  literature  serach  of  the  published  methods  for  calculating 
the  DIDF  revealed  that  the  modified  nonlinear  element  concept,  first  proposed 
for  single  valued  nonlinear  elements  by  Sommerville  and  Atherton  [9],  is 
still  one  of  the  most  widely  used  methods.  It  is  especially  useful  when  the 
secondary  signal  is  a  nonsinusoidal  deterministic  function.  However,  its  use 
in  determining  the  DIDF  for  multivalued  nonlinear  elements  is  much  more 
restricted  than  for  single  valued  nonlinear  elements. 

It  appears  that  the  inverse  DIDF  problem  as  posed  in  Chapter  III  will  be 
difficult  to  apply  in  the  case  of  a  general  nonlinear  element.  This  inverse 
DIDF  problem  consists  of  finding  the  waveform  of  a  secondary  signal  which  will 
result  in  a  specified  modified  nonlinear  element.  This  problem  was  solved  for 
a  restricted  class  of  nonlinear  elements  and  it  was  shown  that  a  general  equa¬ 
tion  for  the  DIDF  of  certain  nonlinear  elements  may  be  found  without  regard  to 
the  waveform  of  the  secondary  signal.  In  particular  it  was  shown  that  only  the 
coefficients  in  the  equation  are  altered  (not  deleted)  when  the  waveform  of  the 
secondary  signal  is  changed.  If  the  order  ot  the  equation  describing  N(e)  is  p 
then  the  order  of  the  equation  describing  the  modified  nonlinear  element  is  not 
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lower  than  p  -  1,  although  it  may  be  larger  than  p.  Thus  when  the  order  of  the 
predominant  terms  in  the  equation  describing  X  (e)  is  large,  the  secondary 
signal  waveform  has  little  effect  on  the  appearance  or  characteristics  of  the 
modified  nonlinear  element.  Such  examples  serve  to  indicate  that  the  use  of 
the  inverse  DIDF  procedure  as  a  synthesis  tool  will  be  of  limited  value.  In  this 
connection  it  was  observed  that  a  very  wide  range  in  the  characteristics  of  the 
modified  nonlinear  element  of  the  relay  function  [N(e)  is  zero  order)  ]  could  be 
affected  by  a  change  in  the  secondary  signal  waveform.  Thus  the  proposal  of 
Oldenburger  and  Ikebe,  in  which  a  relay  function  is  inserted  ahead  of  a  given 
nonlinear  element  in  conjunction  with  the  secondary  signal  of  Equation  (3.20) 
in  this  thesis,  provides  a  versatile  synthesis  tool  with  which  to  change  existing 
nonlinear  characteristics  to  desirable  ones. 

The  modified  nonlinear  element  and  DIDF  curves  given  in  Chapter  IV 
suggest  that  undesirable  nonlinear  characteristics  may  often  be  modified  to 
give  acceptable  characteristics  by  the  proper  choice  of  a  secondary  input  signal. 
These  DIDF  and  modified  nonlinear  element  curves  should  be  beneficial  in 
certain  synthesis  problems. 

The  average  DIDF  method  given  in  Chapter  II  provides  an  alternate 
procedure  for  calculating  the  DIDF.  This  method  reduces  to  a  fairly  simple 
and  compact  form  when  the  nonlinear  element  is  single  valued  and  both 
components  of  the  input  are  sine  waves.  Examples  show  that  in  such  a  case  the 
average  DIDF  method  may  involve  less  computational  work  than  the  modified 
nonlinear  element  method.  The  average  DIDF  method  also  complements  the 
modified  nonlinear  element  method  in  giving  physical  interpretation  to  the  DIDF, 
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although  the  intermediate  step  of  the  latter  method  is  more  useful. 

The  literature  search  revealed  that  very  few  useful  results  have  been 
published  on  effective  procedures  for  calculating  the  DIDF  for  multivalued 
nonlinear  elements.  One  exception  is  the  application  [16,30]  of  the  modified 
nonlinear  element  method  which  is  limited  to  the  case  w'hen  n  >>  1  or  its 
extension  [16]  which  holds  when  n  «  1.  In  this  thesis  the  incorrectness  of 
some  of  the  alternative  procedures  which  have  been  proposed  in  previous  pub¬ 
lications  has  been  demonstrated.  Mahalanabis  and  Oldenburger  [31]  have 
recently  published  a  paper  in  which  they  offer  a  new  procedure  based  upon  an 
"envelope"  description  of  the  dual  input.  In  this  way,  they  obtain  a  single 
expression  for  the  DIDF  which  does  not  include  n  as  a  parameter  but  wliich, 
they  assert,  holds  for  both  n  >  1  and  n  <  1.  It  w-ould  seem  to  this  writer  that 
if  the  modified  nonlinear  element  is  incorrect  for  multivalued  nonlinear 
elements,  then  any  correct  method  must  include  the  frequency  ratio  n  in  the 
analysis  and  final  answer. 

The  average  DIDF  method  offers  such  a  procedure.  It  suffers  from  the 
amount  of  computer  time  required  to  give  a  complete  solution  and  the  number 
of  graphs  required  for  presentation  in  graph  form.  The  TSIDF  of  the  relay 
with  hysteresis  was  calculated  by  this  method  and  presented  in  Chapter  V. 
Unfortunately  the  order  of  integration  cannot  be  reversed  in  the  double¬ 
integral  formulation  [Equations  (2. 12)  and  (2. 13)  ]  when  N(e)  is  multivalued. 
However,  this  is  in  essence  what  is  done  when  the  modified  nonlinear  element 
method  is  incorrectly  applied  to  multivalued  nonlinear  elements  and  results 
in  the  same  answers.  An  approximate  method  for  solving  the  first  or  both 
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integrals  by  analytical  means  would  be  very  helpful.  Further  studies  in  this 
area  might  prove  fruitful. 

Some  experimentally  determined  DIDF's  for  comm  on -multivalued  type 
nonlinear  elements  would  be  of  benefit  to  the  control  engineer.  For  example 
in  pulse  width  modulation  (PWM)  control  systems  incorporating  a  controller 
actuator  which  is  characterized  by  a  relay  with  hysteresis  or  saturation  with 
hysteresis,  it  is  not  uncommon  to  have  an  n  of  6  or  8.  However,  the 
secondary  signal  would  usually  be  a  triangle  wave  in  the  PWM  system. 

The  average  DIDF  method  may  be  extended  to  incorporate  an  input 
consisting  of  the  sum  of  three  or  more  sinusoidal  functions.  Some  useful 
approximate  results  might  be  obtained  in  such  an  investigation  if  one  restricts 
attention  to  single  valued  nonlinear  elements  where  the  parameter  n  may  be 
set  equal  to  unity. 


APPENDIX 


EXAMPLE  CALCULATIONS 
OF  THE 

MODIFIED  NONUNEAR  ELEMENT  AND  DIDF 

The  modified  nonlinear  element  method  as  e:  rplained  in  Chapter  II  is 
used  to  find  the  DIDF  of  the  absquare  when  the  secondary  signal  is  a  triangle 
or  square  wave.  The  average  DIDF  method  developed  in  Chapter  II  is  used  to 
calculate  the  TSIDF.  The  procedures  given  in  the  following  detailed  example 
were  followed  in  deriving  the  DIDF's  tabulated  in  Chapter  IV.  The  cases 
where  the  secondary  signal  is  a  sine  wave  are  the  most  difficult  to  solve  and 
the  average  DIDF  method  provides  some  labor  savings  in  these  cases.  The 
approximate  solutions  for  the  TSIDF  of  dead  band,  limiter,  and  limiter  with 
dead  band  given  in  Chapter  IV  were  found  by  the  approximate  trapezoidal 
integration  rule  in  solving  the  second  integral  of  the  average  DIDF  method. 

A.  Sine  Wave  Secondary  Signal 

The  equation  describing  the  absquare  nonlinear  element  is 

N  <e)  =  e2  sgn  (e)  .  (A-I) 

The  DIDF  is  found  by  Equations  (2.29)  and  (2.31)  which  are  repeated  here  for 
convenience. 
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(A.  2) 


27T 

g(A,B,  p)  =  A  T .  !t  ■  _  I  N(D  sin  cut)sin  cut  dcut 

I )  ta  ^ 


where 


and 


D  £  \Ta2  +  B2  +  2AB  cos  p 


1  27T 

K(A,B)  =  —  fg(A,B,p)dp  . 

Zh  J 

o 

Substituting  Equation  (A.l)  into  Equation  (A. 2)  gives  g(A,B,0)  for  the 
absquare. 


g(A,B,0)  = 


-  B 

1  +  —  cos  p  _  n 

A  2  r  2  •  2 

- ■  —  D  sm  cut  dcut 

77  J 
O 


D 


or 


g(A,  B,  (p) 


8D 

37rA 


(A  +  B  cos  cp) 


The  DIDF  becomes 


4  27r  I - 

K(A,  B)  =  ~ f  (A  +  B  cos  p)  N  A2  +  B2  +  2AB  cos  (p  dp  . 

o 

The  integrand  is  symmetrical  about  p  =  ir  and  is  solved  by  separating  the 
integral  into  two  parts  which  are  easily  identified  as  elleptic  integrals. 


8 


K(A,B)  =  9.  f  n  A2  +  B2  +  2AB  cos  c p  dp 

o7T  A  J 


I 


8  -  * 


B 


3t r2  A 


J  COS  P  n]  1 


A2  +  B2  +  2AB  cos  p  dp 


The  solution  of  the  first  integral  of  Equation  (A.  8)  is  found  to  be 


(A. 3) 


(A. 4) 


(A.  5) 


(A. 6) 


(A.  7) 


(A.  8) 
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Q'=^,A+  B>E(r-k) 


where 


(A. 9) 


2>fAB 
(A  +  B) ‘ 


and  E^k,  y  )  is  the  complete  elliptic  integral  of  the  secoi. 
integral  of  Equation  (A.  8)  is  denoted  by  Q2. 


(A.  10) 


'd.  The  second 


8  B  ? 


Q2  -  37T2  J-  f  cos  <p  J  A2  +  B2  +  2AB  cos  <p  d<p 


(A.  11) 


Q2  =  37r2  sin  $  \l A2  +  B2  +  2AB  cos  <p  )  I* 


sin2  4>  d <f> 


0  \/  A2  +  B2  +  2AB  cos  0 


(A.  12) 


The  solution  of  Equation  (A.  12)  is  shown  by  Byrd  and  Friedman  [38]  to  be 


Q2  i7(A  +  B)(7i?-)[(2-k2)EW-2(k')2F(k)]  ,  (A. 13) 


where 


(k’r  =  1  -  k2 


(A.  14) 


The  modulus  k  is  given  by  Equation  (A.  10)  and  f(-£-,  k)is  the  complete 
elliptic  integral  of  the  first  kind.  Equation  (A.  13)  may  be  written 

Qz  =  97^2~~)  t(A2  +  b2>E  -  (A  -  B)2K]  ,  (A.  15) 


where,  according  to  standard  notation, 


E  A  E(f,  k)  =  E| 


(A.  16) 


^  F(y,  k)  =  F(k)  . 


(A.  17) 


. . 


I 


■ 


Adding  Qt  and  Q2  gives  the  DIDF  of  absquare. 


K(A , B)  =  f  (7A2  +  B2)E  -  (A  -  B)2K] 


(A.  18) 


The  modified  nonlinear  element  N ^Aq,  B^  is  explained  in  Chapter  II.  The 
expression  for  N^Aq,b)  is  given  by  the  integral 


N 


i  27r 

1  r 


(VB)  “J  N(Ao  +  Bsin<!t)#t 


(A.  19) 


Equation  (A.  19)  is  separated  into  two  parts  for  the  absquare;  i.e. , 

2ir-0\ 

S 

tt+Q 


N 


1  7T+01  —  ,, 

(A0>  B)  ~  /  (Ao  +  B  Sin  ^)2d/?t  “27  /  (Ao  +  B  sin  £>t)2d,«t  , 

(A. 20) 


-0 


where 


Ci  =  sin 


_i 


sat 


/A  \ 

I  — I 

\  B  / 


(A. 21) 


N^A^.B^  =-^-  (2Aq2  +  B2)^  +  4AqB  cos  0 }  -  B2  sin  G j  cos  0] 


Evaluating  Equation  (A.  20)  gives 
1 

|^2AQ‘  +  B1-^,  +  4A^B  cos  Gi  -  B‘  sin  0,  cos  G\  I 

J  (A.  22) 

The  DIDF  could  have  been  found  by  using  the  modified  nonlinear  element  of 
Equation  (A.  22)  as  the  first  stage  of  the  two-stage  evaluation  process 
(Chapter  II) .  However,  the  mechanics  of  solving  the  integrals  are  somewhat 
more  involved,  but  the  same  final  expressions  may  be  worked  out.  Equation 


(A. 22)  simplifies  to 

r 


N 


—  M2A  2  +  B2)  sin"1  fl,  +  3A  B 
!  7T  L\  O  '  0 

(vbH 

y  (a2  +yb2)  ;  a  >  b  . 


<  B 


(A. 23) 


a 

is 
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The  input  and  output  signals  are  shown  in  Figure  A.  1.  A  multiplying  factor  c 
may  be  inserted  to  add  to  ihe  generality  of  the  absquare. 


Figure  A-l.  Input  and  Output  Signals  of  Absquare 
for  Computing  Modified  Nonlinear  Element 

B.  Triangle  Wave  Secondary  Signal 

The  application  of  the  modified  nonlinear  element  concept  is  very 
straightforward  when  the  secondary  signal  is  either  a  triangle  wave  or  square 
wave.  In  many  instances  the  "stabilizing"  property  of  these  secondary  signals 
is  not  very  different  than  that  of  a  sine  wave  [32].  however,  the  resulting 
effective  gain  may  be  very  different,  as  was  shown  in  Chapter  III  and 
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elsewhere  [16,33].  The  modified  nonlinear  element  is  first  found  as  follows. 

Assume  the  triangle  wave  crf/T)  and  the  dc  bias  A  shown  in  Figure  A-2. 
The  equation  for  the  triangle  wave  secondary  signal  is 


CT(^t)  =■  B 


■C-*?)  ^  ■ 


-r/2  /°t  £  tt/2 


tt/2  <  ,n  <  —  . 


The  integral  form  for  the  modified  nonlinear  element  is 


H\-b)  -  wf  (A,  -  T-  <’1)  ^ 


(A .  24) 


f  A  -  B 
*  ri,  L  ° 


LJLiiV 

r  /2  / 


(A. 25) 


The  last  integral  of  Equation  (A. 25)  is  best  simplified  su  that 

N- (A„,  B)  =  i /'(a  4/'(a0  (A.  26) 

-flj  NO”'/  N  ' 


x(Aq)  =4  2A^ 


2,  4A°B^2  n  A  8B?.3_ 


(A.  27) 


The  angle  hj  is  found  from  the  relation 


erf''/,)  =  A 


(A.  28) 


»1  “(rH 


(A-29) 


Substituting  for  0 i  in  E  nation  (A.  27)  results  in  the  following  relationships: 
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The  DIDF  is  found  with  the  signal  A  sin  cot  as  the  primary  component  of  the 
input  to  the  modified  nonlinear  element  N  B  ^  ( Figure  A -3) .  It  should  be 
noted  that  Equation  (A. 30)  would  have  to  be  modified  slightly  to  hold  for 
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Figure  A -3.  Modified  Nonlinear  Element  for  Absquare 
with  Triangle  Wave  Secondary  Signal, 

Shown  with  Primary  Input  Signal 


negative  values  of  Aq.  This  will  hold  true  for  most  of  the  modified  nonlinear 
elements  to  be  derived  and  will  be  neglected  hereafter.. 

The  DIDF  is  found  by  using  Equation  (2.5)  which  is  restated  here. 


1  Z7T  ^ 

K(A,  B)  =  —  f  N(  A  sin  cut,  B)  sin  cut  dcut  .  (A. 31) 


Therefore , 


K(A,B) 


4  f  V .  4.  1  A3  sin3  u)t\  . 

~  wA  J  ( AB  sin  +  3 - b - )sm 

o'  ' 


where 


4  r/!'“  /  ;  ,  1  ,\ 

+  tx"  J  \A  sm'  wt  "  T  Bj sm  ^  da;t 


>  i 


K(A,B)  -  2  sin  y,  cos  ■> ,  ^ 


A3/y,  2  sin y ,  cos  y,  sin3  ,  cos  y, 


B  V  2 


4  Id 

A  [4  cosy,  -  — cosJy,)  -  — —  cos 


)  1 


•  i  ,  B  \ 
y,  =  sin"1  sat  (— J  . 


Substituting  this  value  for'/,  into  Equation  (A.  33)  gives 

r 


V  f  A  .  R  \  ^  T 

l] 


t2\  .  ,  /  B\  f  B  /  B2 

B)Sm  [J )  +  J1  -A{3A 


13A 


N  v  / 


A2 

4B 


2B 

A  <  B 


(A. 33) 


(A. 34) 


A  >  B 


(A.  35) 


C.  Square  Wave  Secondary  Signal 


In  this  case  y,  is  always  zero  (Figure  A-l) .  The  square  wave  signal  is 
shown  in  Figure  A -4. 

The  modified  nonlinear  element  is  described  by  the  relation 

*(vB) =  i  Aao  +  B)Id,?l  * SB"  (A°  •  B)  ~h  /  (A„  -  B),dpt  • 

(A. 36) 

Evaluating  these  integrals  shows  the  modified  nonlinear  element  to  be 
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n[ao,b) 


2  BA  ;  A  s  B 
o  o 


A  2  +  B2  ;  A  >  B  . 
o  o 


(A.  37) 
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Figure  A-4.  Bias  and  Square  Wave  Secondary  Signal 


Again  the  DIDF  is  found  by  assuming  the  signal  A  sin  cut  as  the  input  to 
X^A^.b)  and  determining  the  fundamental  output  (Figure  A-3) 


4  P  , 

K(  A,  B)  =—  I  2AB  sin2  cut  dcut 
rA  J 


f  ( A 2  sin2  cut  +  B2)  sin  cut  dcut 


where  again 


i  /  E  \ 
y i  =  sin-  sat/  — 

\  A  I 


(A.  38) 


(A.  39) 


The  DIDF  is  found  to  be 


-it 


K(A,B)  =—  B  sin  sat 

7T 


(!K 


2A  +  YjJ 1  "  sat(f  . 


(A. 40) 


Separating  the  last  equation  into  two  parts  results  in  the  following  simplified 


expressions: 


K(A ,  B)  = 


;  B  <  A 


2B  ;  B  >  A  . 


(A. 41) 


This  shows  the  linear  range  mentioned  in  Chapter  HI  [Equation  (3.31) 
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